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Q ■ Abstract 

. Motivated by the recent interest in dynamical properties of topologically non- 

\ trivial spacetimes, we study linear perturbations of spatially closed Bianchi III 

J> ■ vacuum spacetimes, whose spatial topology is the direct product of a higher genus 

. surface and the circle. We first develop necessary mode functions, vectors, and 

Q \ tensors, and then perform separations of (perturbation) variables. The perturba- 

' tion equations decouple in a way that is similar to but a generalization of those 

of the Regge-Wheeler spherically symmetric case. We further achieve a decou- 
^ \ pling of each set of perturbation equations into gauge-dependent and independent 

parts, by which we obtain wave equations for the gauge- invariant variables. We 
O^' then discuss choices of gauge and stability properties. Details of the compactifi- 

$H . cation of Bianchi III manifolds and spacetimes are presented in an appendix. In 

the other appendices we study scalar field and electromagnetic equations on the 
same background to compare asymptotic properties. 
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^ ! 1 Introduction 

Astronomical observations strongly support the picture of an isotropic and homogeneous 
universe which is expanding at such a rate, relative to its energy density, that it will 
continue to expand forever. Isotropy and homogeneity must be interpreted as valid only 
when the universe is viewed in a coarse-grained sense — on a sufficiently large scale 
that fine details such as clusters of galaxies are blended into a uniform background. For 
this reason most studies of cosmological perturbations have begun with a homogeneous 
and isotropic background and thus with a Friedman-Robertson- Walker (FRW) metric to 
perturb. 

There are however several logical limitations to this point of view. First of all since one 
cannot really observe the entire universe it is perfectly conceivable that the actual universe 
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is only locally but not globally homogeneous and isotropic. For example the negatively 
curved k = —1 FRW models can be spatially compactified in infinitely many ways to yield 
spacetimes which, though indistinguishable from the conventional FRW models on suffi- 
ciently small scales (which might however exceed the observable universe in spatial extent) 
are nevertheless inhomogeneous in the global sense. In particular, at a given instant of 
cosmic time, the length of the shortest closed spacelike geodesic through a given point 
(roughly the shortest distance around the universe from there) would depend upon the 
location of that point in a compactified model. This is a clear (though perhaps unobserv- 
able) violation of global homogeneity. Furthermore since astronomical observations are 
limited to (past) lightlike directions for any given observer, any hypothetical breakdown 
of global homogeneity might be unobservable, even in principle, at a given epoch of cosmic 
time but became observable at a later epoch. 

Both classes of "open" FRW models (the k = —1 hyperbolic and the /c = flat models) 
are amenable to spatial compactification but the A: = models are somewhat unstable 
towards perturbation into locally homogeneous but non-isotropic, Kasner-like behavior in 
which the effective Hubble constant would be directionally dependent and thus (in principle 
at least) observable from any vantage point in the universe. By contrast, the compactified 
k = —1 models are known to be stable in the direction of cosmological expansion and 
thus do not seem to require an excessive "fine-tuning" of the initial conditions (or perhaps 
an inflationary scenario) to protect them from the development of non-isotropic expansion. 

An even more radical departure from the conventional FRW point of view is conceiv- 
able however and even suggested by recent results and conjectures on possible 3-manifold 
topologies and recent insights into the dynamics of expanding cosmological models implied 
by Einstein's equations. Aside from the compact hyperbolic manifolds mentioned above 
which can support k = —1 locally (but not globally) homogeneous and isotropic metrics 
there is a more general class of prime compact 3-manifolds that one can form by "glueing" 
finite volume hyperbolic manifolds across certain special surfaces (called "incompressible 
2-tori") to so-called graph manifolds of infinite fundamental group. This set of graph man- 
ifolds includes for example the trivial and non-trivial circle bundles over compact higher 
genus surfaces. This glueing construction (which is distinct from the connected sum op- 
eration in which the glueing occurs across embedded spheres instead of tori) leads to a 
class of manifolds which no longer admit even an everywhere locally homogeneous and 
isotropic metric since the graph manifolds pieces are incompatible with the extension of 
such a structure from the hyperbolic pieces. 

At first sight it would seem that such manifolds are ruled out, as models for the actual 
universe, by the most basic of astronomical observations. However some recent studies 
of the so-called reduced Hamiltonian for Einstein's equations and their implications for 
cosmological dynamics suggest a different point of view. They support the scenario that, 
under cosmological expansion, the spatial metric (after rescaling by the cube of the mean 
curvature to take out the overall effect of expansion) converges to the hyperbolic metric on 
the hyperbolic components of the manifold but collapses the graph manifold components 
to zero rescaled volume asymptotically. This picture is suggested by the known behavior 
of the reduced Hamiltonian (and its quasi-local variant) which monotonically decreases 
for every non-selfsimilar solution of the field equations and whose value at any instant (in 
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a constant mean curvature or CMC time slicing) is the rescaled volume at that instant. 
The infimum of the reduced Hamiltonian is known to be zero for a pure graph manifold 
of the type we are considering and the asymptotic vanishing of this Hamiltonian corre- 
sponds to the Cheeger-Gromov style collapse of the rescaled volume mentioned above. 
This total "collapse of the graphs" is seen to occur for the explicitly known, spatially 
compactifiable Bianchi vacuum spacetimes of types I, II, III, VIq, and VIII [2] which are 
all graph manifolds of the type we want. By contrast the infimum of the reduced Hamil- 
tonian for compact hyperbolic manifolds is conjectured (through its direct connection to 
the topological invariant known as the cr-constant) to be achieved by convergence to the 
hyperbolic metric and it is known that this canonical metric provides at least an isolated 
local minimum to the reduced Hamiltonian ^ |21 El 13 • 

Thus one is led to consider a model of the universe in which the spatial topology is a 
composite of hyperbolic and graph manifold pieces but for which the Einsteinien dynamics 
predicts that the spatial volume is asymptotically dominated by the (locally homogeneous 
and isotropic) hyperbolic components. Whether a physically plausible cosmological model 
could be developed along these lines (e.g., whether structure formation could be largely 
confined to the hyperbolic components) is not known but, from the mathematical point 
of view, it is worth noting that such manifolds which are composite with respect to the 
torus decomposition described above can still be prime (i.e., indecomposible into more 
fundamental prime manifolds through the sphere decomposition). As such they are still 
among the most basic building blocks of compact 3-manifolds and indeed, taken together 
with the special case of pure compact hyperbolic manifolds, form perhaps the largest 
distinguished subset of such prime manifolds. Together with the (pure) graph manifolds 
of infinite fundamental group they are conjectured to exhaust the family of so-called prime, 
K{Tr, 1) manifolds. If current conjectures (including the still unproven Poincare conjecture) 
on 3-manifold topology are correct the most general compact, connected, orientable 3- 
manifold is constructable as the connect sum of a finite number of K(tt, 1) factors, a finite 
number of x handles and a finite number of spherical space forms (i.e., quotients of 

by finite subgroups of SO (4) which act freely and properly discontinuously on S^). 

Of course there is no known exact solution of Einstein's equations on any of these more 
exotic torus decomposible or sphere decomposible manifolds. However there are known 
solutions available for many of the constituents of such composite manifolds. In particular, 
there are known Bianchi solutions, both vacuum and fluid filled, for most of the graph 
manifold examples mentioned above. Since these could conceivably play an important role 
in the structure and evolution of a plausible universe model we are motivated to study 
their stability properties for much the same reason that people have studied the stability 
properties of the FRW models. 

Going beyond the FRW models and the Bianchi I, V and IX families which include 
such models the simplest manifolds to consider would seem to be those of the Bianchi HI 
(or Thurston's x R) type which can be compactifled to trivial S'^-bundles over higher 
genus surfaces. Starting with a family of explicitly known vacuum solutions we here de- 
velop the tensor harmonics needed for the separation of the spatial coordinates on such 
backgrounds and carry out this separation explicitly for the linearized Einstein equations 
as well as for the Maxwell and massless scalar wave equations. These harmonics have much 
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in common with those developed by Regge and Wheeler ^B] for the study of spherically 
symmetric backgrounds (such as Schwarzschild or Nariai-Kantowski-Sachs solutions) but 
with eigenfunctions of the scalar Laplacian on the higher genus surface playing the role of 
the ordinary spherical harmonics on S'^ in the Regge- Wheeler analysis. For completeness 
though we develop also the harmonics based on the transverse-traceless tensors admitted 
by the higher genus surfaces as well as those based on the (Hodge-) harmonic one-forms 
allowed by such surfaces. Since the first Betti number of 5^ is vanishing and since 5^ 
has a trivial Teichmiiller space these latter classes of tensor harmonics had no analogue 
in the Regge-Wheeler analysis. Furthermore even the scalar harmonics vary with the 
Teichmiiller parameters determining the higher genus surface and thus depend upon the 
particular compactification of chosen in the construction. As we shall see the asymp- 
totic decay properties of the gravitational, electromagnetic or scalar field perturbations 
can vary significantly with the Teichmiiller parameters of the chosen compact surface. 

To obtain a further simplification of the linearized field equations, we also carry out 
a decomposition of the metric perturbations into gauge invariant, gauge dependent and 
constrained variables, the last two subsets being essentially canonically conjugate to each 
other. We then derive the decoupled wave equations satisfied by the various gauge invariant 
perturbation modes and discuss the behavior of their solutions on the chosen backgrounds. 
By choosing a particular gauge for the remaining, unconstrained variables we show how the 
full spacetime metrics of the perturbed models can be expressed in terms of the solutions 
of the gauge invariant wave equations. 

Our analysis lays the groundwork for the study of a number of physically interesting 
questions about these models. For example do the perturbed models decay in some natural, 
invariant sense towards local homogeneity in the asymptotic regime or do the perturbations 
get driven away from asymptotic homogeneity by the Cheeger-Gromov "collapse of graph" 
effect of the background? Does the reduced Hamiltonian for the perturbed solutions 
asymptotically achieve the infimum attained by the background solutions (corresponding 
to the vanishing of the cr-constant for these manifolds) or does it stay bounded away from 
its infimum for a generic perturbations? Can the methods developed here be generalized 
to treat the perturbations of Bianchi II (Thurston Nil), Bianchi VIq (Thurston Sol), and 
Bianchi VIII (Thurston SL(2, R)) backgrounds and, if so, can one answer for those models 
the same questions raised above for the Bianchi III (Thurston x R) models studied 
here? Finally can one recover the essential features of the asymptotic behavior of these 
models from energy arguments for the linearized field equations, which avoid the need for 
an explicit separation of variables and, if so, can one develop these into energy arguments 
for the fully non-linear equations in the case of sufficiently small perturbations? These are 
some of the issues we hope to study in the future. 

The plan of this paper is as follows. In the next section we briefly describe the back- 
ground spacetime we consider and set up the background (canonical) variables. The 
background spacetime is a generalization of the usual Bianchi III model that is spatially 
compactified. Details of the compactification are discussed in Appendix ^ In SecElwe 
construct mode functions, vectors, and tensors on the spatial manifold. Using these in the 
subsequent section we perform separations of (perturbation) variables and write down the 
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(separated) constraint equations. It is claimed that we have four kinds of perturbations. 
The subsequent four sections are devoted to detailed analysis of them, which includes 
decoupling of each set of perturbation equations into gauge-dependent and independent 
parts, and deriving the wave equation for the gauge-invariant variable. In Sects. and 
llfll we discuss choices of gauge. The last section is devoted to a summary with some 
discussions of the stability properties. In Appendices IbI and IHl we study scalar field and 
electromagnetic equations on the same background. All the appendices deal with subjects 
that have not appeared elsewhere and are substantial parts of this paper. 

We employ the abstract index notation j2Sl and use leading Latin letters a, 6, • • • to 
denote abstract indices for vectors and tensors. When useful however we also write a vector 
or tensor without abstract indices, as long as no confusion occurs. Greek indices /i, i^, • • • 
are used to denote spacetime coordinates, running to 3. Middle Latin letters i, j, • • • are 
used for general purposes, e.g. for equations, variables, or generators for infinite groups, 
as well as for spatial coordinates, running 1 to 3. Capital Latin letters /, J, • • • are used 
to label an invariant basis for a Bianchi group, and run from 1 to 3. 

2 The background 

Our spatial manifold M is assumed to be the direct product M ~ Eg x S^, where is 
the closed surface with genus g > 2, and the circle. The geometry is of Bianchi type 
III (or X R in Thurston's terminology). Let M ~ R^ be the universal cover of M. We 
can define on M the invariant 1-forms a (/ = 1 ~ 3) of Bianchi III. Also, we can define 
their duals (invariant vectors) xi""- aXj"" = <^j- In local coordinates they are expressed 
as 

= dx/y, = dy/y, = dz, (1) 

and 

Xi = ydx, X2 = ydy, X3 = dz. (2) 

Since M is compact, the spatial metric Qab on M should be locally rotationally symmetric 
(LRS) to have symmetry (otherwise we cannot make M compact by discrete actions 
of isometry). The background spatial metric is therefore: 

Qab = qihab + q2lab, (the background metric) (3) 

where 

hab = o'^af^b + cr'^a^^'^b (the Standard metric on H"^), (4) 

^ab = c^aC^b (the Standard metric on R), (5) 

and the scale factors qi {i = 1, 2) are constants on M, and will be functions of time t. We 
may view the metric qab as one on the closed manifold M, as well as one on the universal 
cover M. We understand that an appropriate set of identifications of points is specified 
on M when g^b is viewed as a metric on the closed manifold M. Similarly, we may view 
hab as a metric on as well as one on H^, and view lab as a metric on as well as one 
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on R. (For details of the compactification, see Appendix where however we exphcitly 
distinguish metrics on universal covers from ones on compactified manifolds by writing 
tilde on the metrics of the universal covers.) 

We should note that and (or xi and X2) cannot be defined on the compact quo- 
tient M by themselves, since they are not invariant under the isometrics of H^. However, 
they will always appear in a combined form as in the metric (HJ or in the area 2-form 
diifi = (T^ A o"'^, which are both well defined on M. On the other hand, and xs well 
defined on both M and M by themselves. Note that xs is the natural S^-Hher generator 
when viewing M as a 5^ -bundle over T,g. We regard xs = dz not only as a vector but as 
an explicit differential operator acting on the functions on M (or S^), i.e., 

X3:C°°^C°°. (6) 

We also note that in our convention, the scalar curvature for the hyperbolic metric hab 
gives Rh = —2 (in particular, not normalized to give —1). 

We think of qi and q2 as the configuration variables in the background phase space. 
To find the conjugate momentum vr*^^, let fig^ be the determinant of the standard (i.e., 
time-independent) metric qoab = hab + lab with respect to the coordinates {x,y,z). For 
later convenience we also define three-forms dug^ and d^x by 

dfiqo = fJ.qod^X = cr^ A cr^ A cr^. (7) 

The determinant fig of the metric qab is also defined similarly. The two determinants are 
related with fig = (71 y^/U^g . Now, it is easy to see that the conjugate momentum can be 
written in the form 

vr"* = ^igaC^h"^ + 7r2/"*), (the conjugate momentum) (8) 

where /i"'' = Xi^Xi' + X2"X2^ and 1°"^ = Xs^Xs'' are the the inverses of hah and lab, respec- 
tively. In fact, the functions of time tti and 112 are chosen so that they are canonically 
conjugate to qi and q2, as seen from 

©0 = 7^/ d^xTT^^qab = vrigi + 7r2g2, (9) 
J M 

where Cm = / dfig^ is a constant. 

We can calculate the Hamiltonian constraint for the background canonical variables q^ 
and TTj (i = 1, 2), which is 

= /Xg(-^i^2gr^ + l{Tr2fqi^q2 + 2q^^) (10) 
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On the other hand the momentum constraint is found to be triviah 7i"'{q, vr) = —2Dj,7T° 
0. The background Hamiltonian Hq is therefore given by^ 

Hn = -^ I cfxNH 



Cm Jm 

= N^{-Tin^2 + \{i^2?qi\2 + 2), (11) 

where N is the lapse function. With this Hamiltonian we obtain the following Einstein 
equations: 



qi 



TTi 



q2 



d-Ki 

dqi (12) 

-- - N^{'Ki - vr2gf^g2), 
^2 =\Nq2~^'^{2{Tini2 - 2) - Svrig-iga), 

where '= d/dt. We will refer to qi and VTj (i = 1,2) (and N) as the background variables, 
and think of them as given functions of time which satisfy the above evolution equations 
and the constraint = 0. 

The general vacuum solution of Bianchi III is given, up to isometry, by a one parameter 
family of spacetime metrics (§26 of jl7j). for which we write 

9ab = -{t+/t^){dt)a{dt)h + t\hab + {t-/t+)lab, (13) 

where 

= t — k, 

and A: is a real parameter. This metric corresponds to the following solution for the 
background variables: 

qi = tl, q2 = t^t^\ TTi = -2ttl^, 7^2 = -2t+, N = q~^/\ (15) 



^It is well known that one cannot obtain a correct reduced Hamiltonian by the straightforward procedure 
from Bianchi class B models including Bianchi III. This is however not the case for the present LRS case, 
which is spatially compactifiable. (In general, any spatially compactifiable Bianchi models have natural 
Hamiltonians.) 
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3 Harmonic expansion of the perturbation variables 



Since our manifold (M, qab) is naturally a trivial S'^-fiber bundle T,g x any function on M 
may be expanded by the products of mode functions on (S^, hat) and (5^, lab)- This is in 
fact the case when the natural fiber and base are orthogonal to each other. (See Appendix 
^for details of our background solution.) We only consider such backgrounds, and call 
them spatially closed orthogonal Bianchi III backgrounds. ^ In the following therefore we 
first construct each of the sets of mode quantities, i.e., functions, vectors, and symmetric 
tensors, for and for S^, and then construct those for M ~ S"^ x by taking products 
of them. 

Mode Quantities on 

The Laplacian on (S^, lab) is (xa)^ = d'^ jdz^ . The mode functions are therefore, taking 
the period of z-coordinate as 27r, simply given by e*™^^, where the eigenvalue m takes values 

m = 0,±l,--- . (16) 

We refer to m as the fiber eigenvalue. It is convenient to redefine the eigenfunctions 
to be explicitly real so that we can directly calculate perturbation Hamiltonians which 
are second order functions. This can be done by taking the real or imaginary part of 
the complex eigenfunction. Let us define pairs of real eigenfunctions Cm and Cm by the 
relations 

XSCm = -mCm, X3Cm = mCm. (17) 

We assume that these functions are normalized so that the square integrals on {S^,lab) 
give unity: 

/ {Cm)^dz= [ ic^)^dz = l. (18) 

(For example, we can take Cm = {l / \/tt) cos mz , Cm = {1/ ^/^T)smmz.) Mode vectors and 
symmetric tensors can be constructed as Cm,o"^a, Cm^y^a-, Cmlab-, and Cmlab- 
Mode Quantities on Eg 

The mode functions on the hyperbolic surface {T,g,hab) must satisfy 

Af,Sx = -X^Sx, (19) 

where is the Laplacian with respect to hab, and — is its eigenvalue. We will refer to 
A as the base eigenvalue. (For simplicity we omit the subscript A in S\ from now on unless 
it is necessary.) We assume that S is real and normalized so that the square integral on 
(T.g,hab) gives unity, i.e., 

/ dfihS^ = l, (20) 

where dfih = o"^ A cr^ is the area two- form associated with the metric hab- 

^The orthogonal form universal cover metric ^ does not automatically imply this orthogonality of 
{M,qab). See Appendix |X] 
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By analogy with the Regge- Wheeler spherically symmetric case we make "even" 
and "odd" vectors on the from the scalar S when A > 0: 

Sa = —DaS, (the "even" vector on T,g) (21) 
A 

t> = ^Ca'^DbS, (the "odd" vector on S^) (22) 
A 

where Da and eab = (^^[a^'^b] a-^^s, respectively, the covariant derivative operator and volume 
form associated with the metric hab- We always write hats (") on quantities on T,g (except 
hab and eab), and raise and lower indices of hatted quantities (and Cab) by /i"'' and hab, so, 
e.g., ea^ = ^ach^^- The coefficients 1/A are determined so that the square integrals give 
unity, i.e., 

[ dfihh''''SaSb = l, (23) 

and the same for Va- When A = 0, we have Sa = Va = 0, since S = constant. 

Other than these even and odd vectors we need to consider the harmonic vectors Ua- 

ALB^7a = 0, (24) 

where Alb = d6 + 5d is the Laplace-Beltrami operator. This condition is equivalent, in 
the notation of standard differential geometry, to dU = 5U = when the one-form U = Ua 
is on a closed surface, and they can be, in our language using the covariant derivatives, 
expressed as 

D[aUb] = 0, D'^Ua = 0. (25) 

(This is not the same as AhUa = 0.) It is a special case of the Hodge decompositions that 
the space A^{M) of one-forms on a closed Riemannian surface M can be decomposed into 
the direct sum as A^{M) = H^{M) e dA^{M) 5A^{M) with respect to the standard 
L^-norm, where H^{M) is the space of harmonic forms, dA^{M) the space of exact forms, 
and 5A^{M) the space of dual exact forms. It can be easily checked that the spaces of even 
and odd vectors on can be identified with, respectively, dA^{Tjg) and 5A^{Tjg). Hence 
the harmonic vectors definitely complete the possible (smooth) vectors on Sg. Recall that 
harmonic forms are always closed [dtl = 0). From Hodge and de Rham's theorems we 
know that H^{M) is isomorphic to the real- valued cohomology group i/^(M,R). Regge 
and Wheeler did not have to consider harmonic vectors in the spherically symmetric case, 
because the cohomology group H^{S'^, R) of the sphere is trivial. In the present hyperbolic 
case, however, ff^(Sg,R) is nontrivial and given by R'^^, so we have 2g independent 
harmonic vectors on Sg. Let us choose these vectors so that they are L^-orthonormal to 
each other (using the method of Schmidt, if necessary): 

[ dllhh'^'ui'-^UP =6rr', (26) 

where r,r' = 1, • • • , 2g. We call parameter r the cohomology label. All the properties we 
need for the harmonic vectors can be deduced from the definition H25|). which does not 
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depend on r. In this sense r can be viewed as a label that distinguishes "degenerated 
states" . We omit writing these labels when we do not need to specify them explicitly. 

Some of the key relations of the three vectors, which we can easily confirm from the 
definition, are: 



^hYa = -{eig^jYa, eigA^ = < 

and 

b^baXb = -(eig^^) Ya, eig^^ 



A2 + 1, {Ya = Sa) 

A2 + 1, {Ya = Va) (27) 
U, {Ya = lJa) 



A2 + 1, {Ya = Sa) 

1, {Ya = Va) (28) 

1. {Ya = Ua) 



Let us consider symmetric tensors on Eg. They can be split into scalar, vector, and 
tensor parts. The scalar part is the part that is constructed from the scalar mode functions 
S, which part is subdivided into three parts. One of them is the trace part given by Shah-, 
while the others are the (traceless) even and odd parts which are given, when A > 0, by 



Sab = (DabbS + ^Shab ] , (29) 



^-b = \l^^b^aVby (30) 



Again, the coefficients have been chosen so that the square integrals give unity, i.e, 

dflhh''''h'''SabScd = 1, (31) 



and the same for Vab- When A = 0, we define that Sab = ^ab = 0. 

The vector part is the part that cannot be constructed from the mode functions S but 
from the harmonic vectors Ua- It is given by ^ 

Uab = baUb- (32) 

The coefficient (= 1) is as always determined so that the square integrals give unity like 
Sab and Vab- An explicit symmetrization in the right hand side of the above equation is 
not necessary because of the closedness of Ua (See the first equation of Eqs.((2SJ). 

The remaining tensor part is the part that consists of all other possible symmetric 
tensors on {T,g,hab)- From York's decomposition which is analogous to the Hodge 
decomposition, we can identify that they are transverse (i.e., divergenceless) and traceless 
(tt) tensors Wab- 

b'Wab = 0, (33) 
h'^'Wab = 0. (34) 



This tensor is not a zero tensor, as opposed to the flat torus case. 
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{Shab corresponds to York's trace part, while Sab, Vab, and Uab all correspond to York's 
vector part. Beware that York's usage of the term "vector" part (or type) is different 
from ours.) TT tensors appear as variations of the metric hab towards the directions 
that generate Teichmiiller deformations 22 . The freedom of Wab therefore corresponds to 
the linearized Teichmiiller deformations of T,g, which means we have 6(7 — 6 independent 

symmetric TT-tensors on (S^, hab)- We define {s = 1, ■ ■ ■ ,6g — 6) so that they are 

orthonormal to each other, i.e., 

/ dfiHW^l^W^'>' = 6ss'. (35) 

We call parameter s the Teichmiiller label. Like the cohomology label we can think that 
this is a label that distinguishes "degenerated states". We therefore omit writing these 
labels whenever they are not important. 

We can now summarize the symmetric tensors on (S^, hab) as in the table below. It 
is straightforward to check that they are orthogonal to each other with respect to the 
L^-norm. 



X 

Shab 


Even scalar (trace) prt. 


Trace part 


Scalar type 


Generated from 
the eigenfunction 
S = Sx 


Sab 


Even vector prt. 


Traceless part 


Vab 


Odd vector prt. 


Uab 


Harmonic vec. prt. 


Vector type 


Topological 
origin 


Wab 


TT part 


Tensor type 



Table 1: The symmetric tensors on (Tig, hab)- 
Some of the key relations for these tensors are: 



A/iYafe = -(eigA^)^afe, eig^,^ = < 



and 





{Yab 


= Shab) 


A2 + 4, 


[Yab 


= Sab) 


< A2 + 4, 


{Yab 


= Vab) 


4, 


{Yab 


= Uab) 


.2, 


{Yab 


= Wab) 



-3y„6, 

, -2Yab. 

Mode Quantities on M ~ x 5^ 




(36) 



(37) 
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We are now in a position to define the symmetric mode tensors on M. Let S,S £ 
C°°{M) be the products 

S = CmS, (38) 

S = CmS, (39) 

and similarly we make vectors and tensors on M by letting 

Sa = CmSa, (40) 

Va = CmVa, (41) 

Ua = CmUa, (42) 

Sab = CmSab, (43) 

Vab = CmVab, (44) 

Uab = CmUab, (45) 

Wab = CmWab- (46) 

Then, we define our symmetric mode tensors Ei ~ Eg as 

(47) 
(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 

Here, each of Ei to Eq are labeled by the eigenvalues A and m, E^ or Eg is only by m 
with the cohomology label r, and Eg is only by m with the Teichmiiller label s. When 
the Laplacian Aq on M acts on these tensors their response depends in general on the 
eigenvalues A and m, but not on the cohomology labels r and Teichmiiller label s, which 
merely distinguish "degenerated states" for each m, as already remarked. 

For convenience we categorize these basis tensors in two ways, "parts" and "types", 
which we summarize in Table |^ Note that the types are determined with respect to 
the associated quantities on the surface T,g, not to those on the whole M. This usage, 
which seems natural in the present case, does therefore not directly correspond to the one 
adopted in the usual perturbation theory of the isotropic cosmologies. 

For convenience we define a standard L^-inner product in 52 (M) (the space of covariant 
symmetric tensors on M) by 

F,GeS2{M), {F,G)o= [ dfig,q''''q'^FabG,d- (56) 



{Eljab 


=Shab, 




(the even-trace part) 


{E2)ab 


—Saby 




(the even-base part) 


{E'i)ab 


—Slab) 




(the even-fiber part) 


{E/ijab 


=25(aCJ^b), 




(the even-cross part) 


{E^)ab 


=Vab) 




(the odd-base part) 


[Eojab 






(the odd-cross part) 


{E7)ab 


=Uab, 


(the 


: harmonic-base part) 


{Es)ab 




(the 


harmonic-cross part) 


{Eg)ab 


=Wab. 




(the TT part) 
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IT' 


Even part 


Scalar type 


Associated with: 
the eigenfunctions 

on Eg, 
and Cm and Cm on S^, 


771 

-C/2 


TP 


TP 


E5 


Odd part 


Ee 


E7 


Harmonic part 


Vector type 


the harmonic vectors on T,g, 
with Cm and Cm on S^, 


Es 


E9 


TT part 


Tensor type 


and the TT-tensors on T,g, with Cm 



Table 2: The symmetric tensors on (M,qab)- 

The subscript in the bracket is a reminder that the integral is taken with respect to 
the standard volume element HqQ. ^ Contractions are on the other hand taken with the 
time-dependent metric since we stick to the rule that we raise and lower tensors and 
vectors on the spatial manifold M with the metrics q"'^ and Qab- 

The basis tensors Ei defined above are found to be orthogonal with respect to this 
inner product for different kinds: 

{E,,E'^)o = 0, = 1^9), (57) 

where the dash signifies that different values in the labels can be taken. Among the same 
kind we have the following relations: ^ 





'ifiY 




(i 


= 1 ~ 6) 


{E^,Ei)o = < 


ifi)- 


-1 r r 




= 7,8) 






^mm' ^ss' J 




= 9) 



where fi are functions of time given explicitly by 



i 


1 


2 


CO 


4 


5 


6 


7 


8 


9 


h 


2I1 






^qiQ2 













(59) 



4 Perturbation variables and constraints 

Let us apply the basic tools of mode expansion presented in the previous section to the 
gravitational perturbation problem.'' Now that we have the complete set of basis tensors 

^Similarly, we could define another inner product {Ei,Ej), of which integral is taken with respect to 
the time-dependent element fiq. The two are simply related by {Ei,Ej) = qi^/q2{Ei, Ej)o- In this note, 
however, we only use the standard inner product {Ei, Ej)o- 

®For notational simplicity, we assume that A and A' are different even if their numerical values are the 
same if they correspond to distinct eigenfunctions. 

^Before proceeding to the perturbation problem however it is also beneficial to investigate simpler field 
equations on the same background, especially the massless scalar and source-free electromagnetic fields, to 
obtain insights for the properties of the perturbation problem. These materials are given in Appendices IbI 
and 10 
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Ei, we can expand perturbation variables in terms of them. Let 'jab = ^Qab be a pertur- 
bation of the spatial metric qat, and let p"^ = 6tt"'^ be its canonical conjugate. From the 
relations (|58|) we find that we should put 

lab =V7X-^i)a6i 

(60) 

to keep that the components 7* and pi are canonical. In fact, the above relations show 
that 

@^ [ d'xp-'%t, = y^Pif. (61) 

Jm 

Here, 7* = 7*(t) and pi = Pi{t) are functions of time t, and specified by the same labels as 
those of the corresponding basis. The sums are all taken from i = 1 to 9, and all possible 
eigenvalues of and all possible values of other labels (i.e., r and s). 

Constraints 

In general, perturbation variables are constrained by the first variations of the Hamiltonian 
and momentum constraints: 

5n « 0, (62) 
^ 0. (63) 

(See Ref.JJi for a general formula for 571 = DTi.{q, vr)- (7,^) and STi'^ = D'H"'{q, vr)- (7,^) .) 
For our variables they are all obtained by straightforward computations using Eas.()6U|). It 
is, however, worth noting that since the variation of the Hamiltonian constraint itself is a 
scalar (density) only the "scalar type" variables (associated with Ei to E^) can contribute 
to it, and it should become a linear combination of 7' and pi (i = 1 ~ 6). However, the 
odd variables (i.e., 7* and pi for i = 5 and 6) do not actually contribute to this constraint, 
due to the zeros h''^{Ei)ab = l''^{Ei)ab = D''D''{Ei)ab = 0, for i = 5 and 6. Here, Da is the 
covariant derivative operator associated with the spatial (background) metric qab- (See 
Ea. (|229() for a useful formula.) Therefore this constraint takes contributions only from 
the even variables 7* and pi (i = 1 ~ 4). We find after an explicit computation that the 
constraint for each mode is, when A > 0, given by 6TC = fio{qi^,/q2)^^{'DH)S, where 



VH = - I 2m qi + I 7 - Ay — - — 927 




(,.2 IN , 7ri7r2gi 37r22g2\ 3 , ox 4 ^^^^ 

-((A -l)qi + ^ — )7 +2Amgfi7 

- ■T^2qiq2Pl - 92 (TTigi - Vr2g2)P3 ~ 0. 

In the case A = 0, since {E2)ab and (i?4)a6 vanish identically, the above expression is not 
valid, but still the correct expression can be obtained from the same formula by setting 
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7 — 7 — P2 = P4 = 0, as well as setting A = 0. Thus, it is given by 



V 2gi J 



vri7r2gi 37r2^g'2\ 3 (65) 
-1-^1 + ^ 

- '^2qiqiV\ - Q2 (TTigi - 7r292) P3 ~ 0. 

In particular, we denote this function when m = as 

Since the first variation of the momentum constraint is a vector constraint it has four 
kinds of components with respect to the basis covariant vectors Sa, Va, Ua, and Sa^a- Let 
us write the corresponding components as, respectively, VHs, T)Hy, T>Hij, and T>Hi, i.e., 
raising the indices, the variation of the momentum constraint function is of the form 

= -MO [j2i1^Hs)S'' + 5^(Pi7y)F« + Y^iVHuW + E(^^')«2~''^X3"] , (66) 

where the sum is taken over all possible eigenvalues and other labels. By straightforward 
computations the scalar type components, VHs-, 'DHy, and VHi are, when A > 0, found 
to be 



VHs = X{qiPi - 7r27^) - y ^ i'^'il'^ + 2giP2) + m(27r27'^ + ^1^4) ~ 0, (67) 
VHi = -A(7ri7^ + ^2^4) + m(7ri7^ - 7r27^ - 2g2P3) ~ 0, (68) 

VHv = - \l^^{2qiP5 + vra^) + m{qwQ + 27r27^) ~ 0. (69) 

Again, when A = 0, since Ei (i = 2, 4, 5 and 6) vanish identically the above expressions 
are not valid, but still the correct forms are obtained from the same expressions by setting 
the associated variables zeros, at the same time setting A = 0. As a result, VHs and VHy 
are found to vanish identically, while VHi gives a nontrivial contribution: 

^^|A=0) ^ ^^^^^1 _ ^2^3 _ 2^2P3) « 0. (70) 

When m = 0, however, this also vanishes. 

Finally, the vector type component VHu is given by 

VHu = -(2giP7 + 7ri7^) + miqips + 27r27^) ~ 0. (71) 

Note that this is the same as the one obtained by formally setting A = in VHy (and 
changing the indices as 5 — 7 and 6 — 8) . 

The TT-basis Eq does not contribute to any of the constraints, so the associated vari- 
ables 7^ and p9 are not constrained. 
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Decouplings 

Before going into further analysis of the evolutions of the perturbations we present in ad- 
vance how the perturbations decouple. It is direct computations to confirm the following. 

Theorem 1 Perturbations of spatially closed orthogonal Bianchi III vacuum spacetimes 
(with fixed lapse and shift) are excited and evolved independently by either 

(i) the set of the even tensors Ei ~ E4 belonging to the mode for a given pair (A,m), 

(ii) the set of the odd tensors and Eq belonging to the mode for a given pair (A, m), 

(iii) the set of harmonic tensors E'j and Eg belonging to the mode for given m and r, or 

(iv) the TT tensor Eg belonging to the mode for given m and s. 

Here, A is the base eigenvalue, and m is the fiber eigenvalue of the spatial manifold, r 
is the cohomology label, and s is the Teichmuller label. (When two different eigenstates 
have the same numerical value(s) of X and/or m we understand that their eigenvalues are 
distinguished.) 

(When we allow the lapse and shift to perturb, the appropriate kind, i.e., even, odd, har- 
monic, or TT type of mode vectors and scalars are also needed to express the perturbation 
and those kinds of perturbations still decouple from each other.) We call these kinds of 
decoupled perturbations, respectively, the even, odd, harmonic, and tt mode perturba- 
tions. We denote the even mode phase space for given A and m as £\^m- Similarly, OA,m 
is the odd phase space, Um the harmonic one, and Wm the tt one. 

We summarize the mode phase spaces and their constraints in the table below. The 
Hamiltonian for each phase space will be shown in subsequent sections. 



Part 


Mode phase space 


Constraint function(s) 


Even 


Sx,m (A>0,meZ) 


VH, VHs,VHi 


So,m (A = 0,m^O) 




£0,0 (A = 0,m = 0) 


p^(A=m=0) 


Odd 


Ox,m (A>0,mG Z) 


VHv 


Harmonic 


Um (m G Z) 


VHu 


TT 


Wm {m e Z) 






Table 3: Mode phase spaces and their constraints. Each phase space is tacitly assumed 
to be endowed with an appropriate Hamiltonian function, too. These six kinds of mode 
phase spaces should be treated separately. 

For future convenience let us define some terminology here. We call perturbations 
corresponding to all zero eigenvalues (i.e., A = m = for even or odd, or m = for the 
harmonic or tt) zero mode perturbations. We call even and odd perturbations for which 
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either A = or m = 0, or all harmonic and tt perturbations global perturbations. We 
call the generic perturbations that are not zero mode nor global ones local perturbations. 
Perturbations with m = is also called U{1) -symmetric perturbations. 

5 Odd Perturbations 

In this section we deal with the basic problems of the odd perturbations, mainly the 
separation of gauge freedom. The odd perturbations are by definition associated with the 
odd tensors {E^)ah and {EQ)ab- We assume A > as these tensors identically vanish for 
A = 0. In the following we work with each OA,m for a fixed pair of A > and m G Z, and 
assume that the perturbation canonical tensors are given by (See Ea. (|^n)) ) 

where /s and /g are functions of time defined in Eq. (|59j) . The variation of the (covariant) 
shift vector 6Na = S{qabN^) should be an odd vector of the form 

SNa = sVa, (73) 

where s = s{t). We call s the odd shift function. The odd perturbations do not contribute 
to variations of the lapse; 6N = 0. 

The canonical variables which span the odd phase space Oa,™ are given by 7^, 7^, ps, 
and pq. As we have seen in the previous section the odd variables are constrained by a 
single constraint Tiy ~ 0, which defines a three-dimensional subspace in CA,m- Moreover, 
this constraint function generates gauge transformations in this subspace, so the true (i.e., 
gauge-invariant) dynamics is lying in the two dimensional space which is obtained by 
contracting those gauge orbits. We denote this two dimensional phase space as OA,m, and 
call it the gauge-invariant odd phase space. We call the canonical variables in Oa,™ the 
gauge-invariant (canonical) variables, and denote them as Q and P. ^ 

One of the main tasks in this section is to find a canonical transformation such that 
it defines the gauge-invariant canonical variables, and then apply it to the original Hamil- 
tonian in Oa,™,, which is presented right below. As a result we obtain a gauge-decoupled 
form of Hamiltonian -ffoddi for which the gauge-invariant part H^J^^{Q,P) determine the 
dynamics of the gauge-invariant variables. We present the wave equation for the gauge- 
invariant variable Q. The evolution of the other gauge variable will be solved in terms of 
the gauge-invariant variables in SecEl 

The evolution equations and Hamiltonian for the odd variables 

The evolution equations for the odd variables can be reduced from the general equations 
of motion by substituting Eqs. (|7^ and (|7!^ . See Ref.jE] for the general perturbation 

^Although a constraint function is also gauge-invariant, we will, unless otherwise stated, only refer to 
Q or P as a gauge- invariant variable. 
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Hamilton equations in the case of = 1, Na = 0, and 6N = = SNa- Note that, since 
the lapse function N for our background is not unity we need to generahze the formula by 
putting d/dt {l/N)d/dt. We also add the contributions of the perturbed shift. After a 
straightforward computation we obtain: 



7 



P5 



Pa 



N 



N 



+ V2(A2 + 2)s, 



N 



7r2g27 +qiq2P& 
1 



2 , "1^ + VTi + -vrlg^ - ■KiTT2qi ^q2 ) 7^ 



— m 



A2 + 2 



-7 



(vTigi - ■K2q2)PZ 



A2 + 2 _i 



m 



A2 + 2 5 



(\2 -^2— 
A + 7ri7r2 + -712^1 q2 I 7 



T^2q2PG 



+ 2rmT2qi s. 



(74) 



We can easily reconstruct the Hamiltonian -ffodd(7^)7^)P5)P6) for the odd variables from 
the above equations of motion; It is ^: 



N 



odd 



qiwq^ 



1/2 2 ^2 —2 2 —1 \ / 5\ 

- ( m + VTi + -712^1 92 - 7ri7r2gi 92 I (7 ) 



/A2 + 2 5 6 

+ m\i ^—Yr 
+ qipl + \q1q2pl 

+ (vTiQi - ■K2q2)l^Pb + T^2q2l^PQ 



-VHv, 

qi 



(75) 



where T>Hv = VHyij^ ,P5,Pg) is the constraint function for the odd perturbation l|69|) . 

We comment that the terms in the shift function s in the Hamiltonian and evolution 
equations can easily be added after computing them with s = 0. The easiest way to do 



®One can in principle get the same Hamiltonian by computing the second variation of the Hamiltonian 
constraint function. (See, e.g., Ref . [141 for a general formula of the second variation.) It would, however, 
require tremendous calculations, due to the second order terms. 
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this is to note the fact that the perturbation Hamiltonian is obtained by computing (half 
of) the second variation of the background Hamiltonian: 



H = -65 Ho 
2 



= ]-66 [ d\{NH + NJi'^) (76) 

= [ d^xi^NSdn + dNdn + SNadn"). 

JM 2 

(We dropped the term NaSSTC^, which is zero for our background with Na = 0.) The first 
term in the integration contributes to the Hamiltonian for which s = 0. The second term 
is zero for the odd perturbations (but not for the even ones). The last term is the one we 
want to calculate, which is 

/ d^xdNaSH" = - [ duosVaiVHv)^ 

JM JM 

= -sq^'VHv [ dMcn^fh^'VaVb C^^) 
Jm 

= -sq^^VHy. 

Thus, we obtain the final form of the Hamiltonian with an arbitrary s. (Note that we 
already have the explicit form of T)Hy.) The evolution equations with s may be obtained 
from that Hamiltonian. This technique will be more helpful for the even perturbation case 
where much longer computations are required. 



Canonical Transformation 

To decouple the Hamiltonian into gauge-dependent and independent parts and obtain 
the gauge-invariant phase subspace Oa,™, we need to find a canonical transformation such 
that one of the new momentum variables coincides with the constraint; let (7*,p*) be 
the canonical pair thereof, for which the momentum is constrained = 0. Once such a 
canonical transformation has been found, the other canonical pair, denoted as (Q,P), au- 
tomatically becomes gauge-invariant jl2lll4j . We find such a canonical transformation by 
the method of generating function (also JH)- Let 5(7*, Pi) be the generating function 
we are looking for, which is a function of the original configuration variables 7' = (7^, 7^) 
and new momenta Pi = {P,p^:). Since both the original and transformed perturbation 
equations should be linear, it is sufficient for us to consider linear transformations jllj . 
We can therefore assume that S is second order: 

i=5 j=5 1=5 

where Uij = Ci{ij)- Each (3i should be a linear function of Pj, while aij depend only on the 
background variables. 
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We remark that we can redefine the constraint function as 

Cy = —VHv, (79) 

using a (nonzero) arbitrary function of time 9 = 9{t). We call this parameter function the 
mock shift function. The factor qi in the above definition is just for convenience. This 
redefinition is in fact possible, since the original constraint surface specified by T>Hv = 
is equivalent to that for Cy = 0. Since we define a new momentum variable p^, to 
coincide with Cy, this function is one of those parameters which control the canonical 
transformation we consider. 

We set the following Hamilton-Jacobi-like equation 

dS 

p. = Cy{^\pi = -Q-). (80) 

Substituting Eq. (|78|) . we obtain the following equations for the generating function coef- 
ficients aij and /Jj: 



e(V2(A2 + 2)/35-m/36)=P*, (81) 



V2(A2 + 2)a56 + m{am + 2T:2qi^) = 0, (82) 

— ^ — (2a55 + 7rigi )+ma56 = 0. (83) 

These equations do not completely determine aij and /3j; for example, we can regard 
(3q and aee as free functions. This freedom corresponds to the freedom of performing 
canonical transformations among the gauge-invariant variables like {Q,P) — > {Q',P'). 
As we will see below, (3q solely determines the configuration variable Q, and once Pq is 
fixed, then age becomes a parameter that controls the canonical transformations of the 
form {Q,P) — > {Q,P'). We exploit this freedom to make the final form of the gauge- 
invariant Hamiltonian neat, especially "diagonal". For later convenience we "solve" the 
above equations for Oij in terms of a single oqq: 

, -In 1 -1 



056 = "^("66 + 27r2Q'i 



where we have defined a normalized fiber eigenvalue, 

m 



(84) 



(85) 



We simply set 
Then from Ea. (|81() . we have 



\/FT2' 

f36 = P + p,. (86) 




20 



The gauge-invariant variable Q, conjugate to the momentum P is now obtained from the 
generating function as 

^ dP ^ dP^ 

On the other hand, the gauge variable 7=,,, conjugate to is given by 



- ^ - ^ i 
dp* dp*: 

1=5 

( V 1 



(89) 



^/2 V2(A2 + 2)e 



7^ + 7^ 



The inverse is therefore 

,,5 



7" = \/2(A2 + 2)0(Q-7,), 
76 = (1 -m0)Q + m07,. 



(90) 



We also obtain the following transformation between the original and new momentum 
variables: 

dS 5 6 ^ 

Ph = ^ = "557 + "567 + P5 

U 



A2 + 2 ( V 1 \ (91) 



+ V 7, T^l(ll ^^7* + 

dS 5 & n 

P6 = = "567 + "667 + P6 



^ V2(A2 + 2)0 



= (q66 + 2rm:2qi ^0)Q + P - 2rmT2qi ^^7* +P*- 

Equations (|91j) and (|9()|) complete the canonical transformation that takes the Hamiltonian 
into a desired gauge-decoupled form. The functions 9 and age are free parameter functions 
that can be specified at any time one wants. 

It is easy to see that the variable Q = (yl\pi)y' +7^ is certainly gauge-invariant. See 
the induced map (|169|) . 

The Gauge-Decoupled Form of Hamiltonian 

The total Hamiltonian //odd for the new set of canonical variables should be of the following 
(semi-) decoupled form 

H,^d = H^^^{Q,P) + Jp*, (92) 
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where J = J{Q, P, 7*,p*, s; 9) is a function that possibly depends on all the new variables 
Q, P, 7*, and p*, and the shift function s (and also parametrically on the mock shift 9). 
We call J the multiplier function or renormalized shift function for the odd variables. This 
form of Hamiltonian is a consequence of our system being a first-constrained system (e.g., 
121 )• In fact, the Hamilton equations for the gauge variables 

7* = J + ^— p*, 

P* = -t;— P* 

97* 

are consistent with the constraint 

p, = 0. (94) 

Moreover, we can easily see that the evolution of the gauge-invariant variables Q and P is 
determined, up to constraint, only from the gauge-invariant part Hamiltonian H^^^{Q, P). 
In fact, while we have 

^ dP dP ^ dp^*' 

dQ dQ dQ^* 

from the total Hamiltonian, the last term for each equation vanishes if the constraint 
= is imposed. 

We can find the decoupled Hamiltonian -f^odd(Q) -fj 7*;P*) from the original Hamilto- 
nian Hoddil^ ,1^ ,P5tP6) using the definition of the new variables (|55)) . (|S^ and Since 
our canonical transformation is time-dependent, however, we must add the time-derivative 
of the generating function S according to the standard prescription. The correct relation 
is therefore 

dS 

Hodd = Hodd + -Q^, (96) 

where the time derivative operator {d/dt) acts only on the background variables (i.e., they 
do not act on the perturbation variables like Q or P). Prom a direct computation we find 
that the the gauge-invariant part H^J^^{Q,P) is given by 

+ 2TT2{2u'^qi + q2)am + qiu{a^&f\ IT^T^ (97) 

1 ^2 NuP"^ NQP 

+ o"66<y + o , — + "^odd 1=, 

2 2^ qiy/q2 

where we have defined a combined scale function 

u = u'^qi + q2, (98) 
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and a cross term function 



Codd = (2z^ qi + ^2) + Qiu aee, (99) 

for convenience. We have used Eqs. (|Hl|) to express aij in terms of age only, so age is 
now considered as a free function of time. Also, observe that the mock shift and shift 
function s do not appear in H^J^^ as they should not. 
As for the multiplier function J, we find that 

JiQ,P,l*,P*,s;9) = i{2TT2 +91066) (^u- ^2"^ 2 'j) qu/q^ ^) ^ 



/ m qi\ NP 

+ 7T 7* 



(100) 



+ {u + 



Qi 

qi{l - 2me)\ Np^ s 



(A2 + 2)6*2 ; 2^ 9' 



This function determines the evolution of the gauge variable 7=,, through the equation of 
motion 

% = J\o = JiQ,Pn*,P* = 0,s;9). (101) 

Note that we can think of the gauge-invariant variables Q and P in the above equation 
as given functions, since the evolution of the gauge-invariant variables is determined in- 
dependently of the gauge variable 7=,, . We will need the above equation to determine the 
actual evolution of the metric functions 7^ and 7^, especially the evolution expressed in 
terms of Q and P (or Q). 

It is interesting to note that even if we leave the shift function unperturbed (i.e., s = 0) 
we can arbitrarily specify the "gauge velocity" J by suitably specifying 9{t), and thereby 
we can virtually have any desired profile of the gauge variable 7*(t). In this sense we 
can say that the mock shift function 9 mimics the role of the shift function s. We will 
see however that this feature is superficial in the sense that the profiles of the metric 
functions 7^ and 7^ are not affected by the choice of 9 as long as s = 0. In other words, 
this arbitrariness is an artifact contained in the canonical transformation we consider. If 
we allow the shift function s to vary, however, the gauge function 9 automatically gets 
related, after solving the gauge equation of motion, to s. (cf. Sec. El) 

The Wave Equation for the Odd Perturbation 

To eliminate unnecessary ambiguities, let us choose 

a66 = ^ - 2) (102) 
qi \ u J 
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to make Codd = 0, with which the gauge-invariant Hamiltonian H^J^^ becomes a diagonal 
form. In this case, the equations for the gauge-invariant odd variables (up to constraint) 
become 



■ NuP 



VQ2 

u 4uqi u^qi J qi^ 

Then, eliminating P we obtain the following wave equation for Q: 

N _ f -TTiqi + 3TT2q2 , VTiQ'i - 27r2g2\\ A 



(103) 



Q _ A^v . 

\N V 2^192 uqi J 



+ ^ f (A2 + 2)u -q2 + ^(5^292 - 27rigi) + ^(vngi - 27r2<72)l Q = 0. 
9192 V uqi J 



(104) 



Here, we have used the background Einstein equations (fT^ to express the time-derivatives 
of the background variables. This formal equation is valid for any spatially closed orthog- 
onal Bianchi III vacuum background (with vanishing shift). (We also remark that since 
the definition of Q is independent of the choice of aij, the above wave equation for Q does 
not depend upon the choice 1)102(1 .) 

If we substitute the exact solution ()15j) into the above wave equation we obtain an 
explicit form of wave equation: 



where u = z^^ti -|- ^ 



6 Harmonic Perturbations 

We see in this section that the harmonic perturbations can be (formally) viewed as a limit 
of the odd perturbations and because of that, we are required to perform almost no new 
calculations on top of those for the odd ones. 

Remember that the harmonic perturbations are those associated with the harmonic 
vectors on Tig. The corresponding basis tensors are Ej and i^g, which are specified by 
the fiber eigenvalue m. (More precisely, for each m there are 2g independent sets of E-j 
and £"8, in accordance with the 2g independent harmonic vectors C/i^^ (r = I,-- - ,2g).) 
In the following we work with each mode phase space Um = l^m^ for a given m (and r). 
(We hereafter omit writing the cohomology label r for simplicity). We assume that the 
canonical perturbation tensors are given by (cf. Ea. (|^n]) ') 

lab=lHEr)ab + lHEs)ab, (106) 

P''' = l^,o{f7P7{E7r' + fsPsiEsr"), (107) 
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where the functions /y and fs are defined in Eg . (1591) . Observe that Um is spanned by 7'^, 
7^, pj, and ps, and therefore dim^m, = 4. The variation of the shift vector 5Na should be 
a harmonic type vector of the form 



where s = s{t). We call s the harmonic shift function. (Using the same notation of the odd 
shift function s may not cause confusion.) The harmonic perturbations do not contribute 
to variations of the lapse; 6N = 0. 

The variables are constrained by a single constraint T)Hu (See Eq. (|7H) ). A similar 
count to the odd case gives 4— 1 — 1 = 2 as the dimension of the gauge invariant phase 
space Um- 

What we have to do to obtain a gauge-decoupled Hamiltonian and thereby obtain a set 
of gauge invariant equations is completely parallel to those for the odd perturbations. Re- 
member that both the odd symmetric tensor and harmonic symmetric tensor are traceless 
(on Eg); h°'^Vah = h"'^Uab = 0. Notice also that almost all relations for the harmonic vector 
U"" related to the covariant derivative are the same for the odd vector V"" if we assume 
A = 0. For example (see Eq.^), AhU" = -U" is the same as AhV" = -(A^ + l)^'' if we 
put A = 0. (We are talking about "formal" resemblances here, forgetting the fact that the 
odd vector is not defined for A = 0.) The only difference is that the odd vector does not 
satisfy the closed condition that the harmonic vector does. That is, we have D\^gU^ = 
(see Eq.(j2ni)), but D^gV^ / 0. This difference however does not affect our calculations in 
any way, since these vectors contribute to the perturbations only through the symmetrized 
tensors Vab oc -D(aV'b) or tJab = Di^JJy^{= DaUb)- Due to these properties, all the necessary 
formulas (e.g., Hamiltonians) for the harmonic perturbations can be obtained from the 
ones for the odd perturbations by (formally) putting A = 0, as well as changing the indices 
5^7 and 6^8, if necessary. So, we are not trying to re-present most of the results. 

One of our interests may however be the wave equation for the background solution 
(|T5|l . which is given by 



where u = {m? /2)t]^ + (This can be obtained from Ea. ()lU5() by putting A = and 

therefore v = m/y/2.) Here, the gauge-invariant variable Q for the harmonic perturbation 
is defined by (cf. Ea. (|HH|) l 



5Na = sU, 



(108) 




(109) 




(110) 



When m = (the zero mode), this equation becomes 



for which the general solution is given by 



^ + 72^ '3 = 0, 



(111) 




(112) 



where Ci and C2 are integration constants. 



25 



7 Transverse- Traceless Perturbations 



Remember that the variables for the transverse-traceless (tt) perturbations are not con- 
strained, so we do not need any extra steps to get gauge decouphng. 
The TT-perturbation variables can be written 

lab = l^iE9)ab, 

We omit writing the subscripts s that distinguish the 6g — 6 independent TT-tensors. We 
remark again that each mode (for a fixed s) is parameterized by the fiber eigenvalue m, 
but not by the base one A. The tt perturbations do not contribute to perturbations of 
the lapse function and shift vector; 6N = 0, 6Na = 0. 

It is not difficult to find by a straightforward calculation the evolution equations for 
the TT-variables 7^ and pf, which are given by 

7' = A^(2(Z?^^ + (7rigi-7r2<?2) ^ 



AT ( / 2 , 2 , 1 2 -2 2 -I \ 1 I ^ \ P9 

pg = -N (m + + -TT2qi q2 - vri7r2gi q2)- — + (vrigi - 7T2q2)- 



2 2gi^ qi^ 

Let us define for coherent notation 

Q = ^9, P=P9. (115) 

Eliminating P, we find that the wave equation for Q is given by 



Q- I ^ + „ ^^ (vTigi - 3TT2q2) I Q 




1 5 

'm^q2^ - + 2^i^2 - -^-^211^12) ) Q = 0. 



(116) 



Here, we used the background Einstein equation 1)12(1 to eliminate time-derivatives of the 
background variables. 

Note however that since 7^ contributes to the perturbation of metric in the form 
^^ab = l'^{Eg)ah and the base {Eg)ab is tangent to the base ((-£^9)06X3'' = 0)i we can say 
that 7^ belongs to the 'base' (not 'fiber') part, and the base part of the background metric 
has unique scale function qi. A more natural choice of variable is therefore to take the 
rescaled variable^*^ 

Qb^iV<1i = QI<1i- (117) 
(The character B stands for 'base.') Let us rewrite the wave equation for Qb; which is 
given by 



Qb- + ^ ^(vri(?i + TTags) Qb + N\m'q2')QB = 0. (118) 




^''Similar arguments for the other kinds of perturbations are not trivial, since a gauge-invariant variable 
is in general a hnear combination of variables that belong to different parts 'base', 'fiber', and 'cross' ones. 
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It is interesting to observe that this wave equation for the TT perturbation is exactly the 
same as the wave equation for the scalar field in case of X = 0. (See Eg . (1242(1 in Appendix 

El) 

On the exact background solution (|15() . we have 



2i ■ o ^4 



Qb + -=—Qb + ^^Qb = 0. (119) 



Solutions in some special cases 

(i) The flat background {k = 0) solution with m 7^ is 



Qb = - {Ci sin mt + C2 cos mt) , (120) 



where Ci and C2 are integration constants. 

(ii) The zero mode (m = 0) solution with arbitrary k is 



fCi + Calog^, (k^O) , ^ 

\Ci + C2\. {k = 0) 



8 Even Perturbations 



The remaining perturbations are the even perturbations, which are generated by the even 
mode tensors {Ei)ab (i = 1 ~ 4). It is worth bearing in mind that, although any kind of 
computation becomes much longer and sometimes hideous due to the increased number of 
basic variables and the complicated couplings among them, the basic procedure to get a 
decoupled Hamiltonian and wave equation is completely systematic and the same as that 
for the odd case. (We however do not introduce mock shift and lapse functions similar to 
the odd mock shift function (cf. Ea. (|79() ) to simplify the computations and arguments.) 

Remember that two of the basis tensors, {E2)ab and {E4)ab, identically vanish when 
A = 0, so this case should be treated separately. We say that the perturbation is generic 
when A > 0, and nongeneric when A = 0. We first deal with the generic case. 

8.1 Generic (A > 0) perturbations 

In the following we work with the generic even mode phase space <?A,m for given A > and 
m G Z, and suppose that the perturbation canonical tensors are given by (see Eq. (|6flj) 'l 

4 

lab = ^7'{Ei)ab, 

(122) 

Pab = f^qo^fiPi{EiT^, 
i=l 
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where the functions fi (z = 1 ~ 4) are defined in Eg . (1591) . The eight (canonical) variables 
7*, Pi (i = 1 ~ 4) span the mode phase space Sx^m, so dimSx^rn = 8. The variation of the 
lapse function and shift vector should be of the form 

6N = soS, 

- , (123) 

5Na = SiSa + SsSa a, 

where Si {i = 0, 1,3) are functions of time. We call sq the even lapse function, si the even 
base shift function, S3 the even fiber shift function. The gauge-invariant dynamics is of 
course not affected by these functions. 

As we have seen in Sec El the even variables are constrained by three constraints 
DH ~ 0, T)Hs ~ 0, and VHi ~ 0, which define a five (= 8 — 3) dimensional subspace 
in £x,m- Moreover, each function generates independent gauge transformations in this 
subspace. The gauge-invariant dynamics is therefore in a two (= 5 — 3) dimensional 
subspace, which we denote Ex,mi and call the gauge-invariant (even) phase space. We call 
the canonical variables in £x,m the gauge-invariant (canonical) variables, and denote them 
as Q and P. 

Before performing a canonical transformation to get a decoupled Hamiltonian, we need 
to find the evolution equations and Hamiltonian for the original variables. 

The evolution equations and Hamiltonian for the even variables 

It is a lengthy but straightforward computation to obtain evolution equations, using 
Eqs. (|122j) . For simplicity we only present the evolution equations for the vanishing even 
lapse and shift functions (i.e., synchronous gauge) sq = si = S3 = 0. (The general equa- 
tions with arbitrary even lapse and shift functions may be easily recovered from the general 
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Hamiltonian given after the following evolution equations.) 



N 



3 



Q1VQ2.1 1 

^ 7 = -2^2517 - Qiq2P3, 
^ 7 = (TTigi - vr2Q2)7 + 2giP2, 

91^92.3 -1 2 1 1^ Q \ 3 

^ 7 = -T^2qx I2I - i^y^m - 3vr2g2)7 

- 92(^1^1 - q2P'i), 

^ 7 = 7r2g27 + Q'l92P4, 

. (2-2 2 2\ 1 , 1 /'r,\2 ,0 2-1 \ 3 

^ Pi = -(7r2'7i 92 - "T- )7 + 4 (2-^ + 3vr2 ^2) 7 

- mA7'^ + vr2gf ^g2^P3, 
Q'i\/^ ■_l/^2, 2 -1 ,12-22 



^ V2 = -i^\rn ^T^X -T^\T^2qx 92+2^29l 92)7 



;A2 + 2 /A 3 4 

+ \/ — 7, — 7r7 ~ "^7 



2 V2 

- (TTigi - ■K2q2)P2, 



. 1 /„.2 , o 2 -1 A 1 , ^ /A2 + 2 2 

1/^3 2 1 -1 -A 3 

-2(4^2 + 2^1^291^2 -9192 1 7 

1 1/ 

+ 2^29iPi + 2 (^1 91 - 3vr292)p3, 



91^/9^. , 1 / A2 + 2 2 
-p4 = — mA7 — my 7 



^7r29i ^ (27ri9i + vr292)7'' - 't^2 92P4- 



(124) 



It is at once possible to read off the Hamiltonian which gives rise to these equations, 
and also easy to generalize it to include the even lapse and shift functions (see the comment 
below Eq.((7SI)). The Hamiltonian with arbitrary even lapse and shift functions is given 
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by 



N 



Q1VQ2 



1/ 2 -2 2 2\ ^ 1\2 

o {^2Ql 92 - "1 ) (7 ) 



7ri7r2gi q2 + ^'^2^1 92 ) 17 ) 



2\2 



1 / 3 2 1 -1 
+ 4 i 4^2 + 2 ^1^29192 

+ (2vrig'i + 7r2q'2) (7'')^ 

- ^ (2A2 + 37r22g-ig2) 7S' 



9192 ^ 



3\2 



A2 + 2/Aq 4\ 2 

7 — m7 7 



, 22, 22 
+ 9iP2 + 92P3 

+ ^9192 {pI - 2piP3) 

- '^2qi^qWp3 + (TTigi - TT2q2h^P2 

- ^T^2qil^Pl - ^(TTigi - 37r2q'2)7^P3 + T^2q2j^P4: 

-J^VH - -VHs - -VHi, 
9i\/92 91 q2 



(125) 



+ 



where DH, T>Hs, and "DHi are the constraint functions for the even variables, given in 
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Canonical Transformation 

We need to find a canonical transformation so that three of the new momenta coincide 
with the three constraint functions. The remaining new momentum is identified with 
the gauge-invariant momentum P that is conjugate to the gauge-invariant configuration 
variable Q. As in the odd case we write the generating function in the form 

1 ^ ^ 

= 2 E«^^•^V + E/5^^^ (126) 

i,j=l i=l 

where Oij = ct{ij)- The coefficients Oij are functions of background variables only, while 
Pi depend linearly upon P. For convenience, we redefine the constraints in the form 

We remark that the constraint functions in hnear perturbation theory automatically strongly commute, 
which is the condition that is necessary (and sufficient) to assure that they can be new momentum variables 
after a canonical transformation. In fact, since the constraint functions are linear, their brackets do not 
depend on the perturbation variables, and [3 since the brackets at least weakly (i.e., up to the constraints) 
commute with each other, this in turn implies that they must strongly commute. 
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Pi = Pi — Ki {), (i = 1,2, 4), where Ki is a function of 7^, 7^, 7^, 7^, and p^. (This 
is actually possible.) Although this redefinition is not necessary, it helps us to make 
computation more transparent and easier. Let 

H} = {llllQ,lt) (127) 
be new configuration variables, and let 

{Pi*} = {pu,P2*,P,P4*) (128) 



be their conjugate momenta. The equations coming from the Hamilton-Jacobi-like equa- 
tions, PiiY,Pi = = Pu, P2{.l\Pi = = P2*, and P4i-f\pi = = Pi* are the 



following: 

(i) From Pi = cio + cii7^ H h 0147"^ = pu: 

cio = /?i + Ai/33=pi*, (129) 
cii = an + Aiai3 + + ^ + ^ = 0, (130) 

T^2qi T^2q2 2gfi^ 

C12 = ai2 + Aiq;23 + ^ /^^^ = 0> (131) 
V2 7r2 qi 

- -uA , 2 (2(A^-l)+7ri7r2) gi-37r2^g2 „ 

ci3 = aia + Aiaas H ^ = 0, (132) 

47r2gig2 

2 A TO 

ci4 = ai4 + Aia34 = 0. (133) 

7r2g2 

(ii) Prom P2 = C20 + C2i7^ H h 0247^ = P2*- 

C20 = /?2 + A2/33 =P2*, (134) 

, -2TO^7ri7r2gi^ + A^(2A^gig2 + Wgi^ + 7r2^g2^) „ 

C21 = ai2 + A2ai3 + ^ = 0, (135) 

2^2(2 + A2)A7r2gi2g2 

- , A , ^^ + 7ri7r2 , . 

C22 = "22 + A2Q!23 H 7: = 0, (136) 

27r2gi 

, 4TO2 7r2'gi + A=^(2(2(A^ - 1) + ^1^2)^1 + TTaV) 

C23 = 023 + A2a33 + /o/o ^ \2A X = ^' ^^^^^ 

4Y^2(2 + A^)A7r2 gi <?2 

TO (-2 A^ gi + 7r2 (^1 gi - 2 7r2 92)) . . 

C24 = 0:24 + A2a:34 H ^ , , = 0. (138) 

V2(2 + A2)7r2gi(?2 
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(iii) Prom P4 = C40 + 0417-*^ H h 0447^ = p4*: 



_ 2mp3 

C40 = P4-\ — = P4*, (139) 

2mai3 "iTTi 

C41 = ai4 H ^ = 0, (140) 

2ma23 ^ 

C32 = a24 + = 0, (141) 

2 m 033 m7r2 

C43 ^ a34 H h -; = 0, (142) 

A A 52 

_ 2mQ;34 TTi 

C44 = 0:44 -I r \ = 0. (143) 

A q2 



For convenience we have defined 

A,.!Ii-^ = S-2!^, 
7r2 qi A^ 

A2 = ^ V, 
V2(A2 + 2) 

where 

7r2 gi A2 



(144) 



S.^-^ + 2!:V- (145) 



Let us first look at the set of equations Cjo = Pi* (i = 1, 2, and 4), wliicli consists of three 
equations for four unknowns /3i (i = 1 ~ 4). To determine them we put = ps* = P, 
where P is the gauge-invariant momentum, then we immediately obtain the following 
solution 

(3i = pu - AiP, 

02 = P2* - A2P, 

Ps = P, (146) 
2mP 



Pa = PA* 



A 



The remaining 12 equations, cu = C2i = c^i = (i = 1 4), for 10 unknowns 
Oiij = a(^ij) = 1 ~ 4) look overdetermined at first, but this system is underdetermined. 
In fact, given, say, 033, we find we can solve in terms of 033 all equations except three 
equations C21 = C24 = C41 = 0, and these three equations are automatically satisfied as long 
as the background Einstein equations are satisfied. So, the equations are an underdeter- 
mined system. In later calculations we will express Uij in terms of 033 when it is useful. In 
such a case, 033 is regarded as an arbitrary (prescribed) function oft. ("Prescribed" means 
that a33 should be prescribed before solving the perturbation equations.) As pointed out 
in the section for the odd perturbations, this is a result of the fact that there remains 
freedom to canonically transform among a given gauge invariant pair of variables. 
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We can determine the gauge-invariant variable Q from the generating function S with 
the sokition (|146|) . We get 

Q = - 
^ dP 

= -Ai7^-A27' + 7=^-^^. (147) 

A 

Let us denote the remaining gauge configuration variables as 7* = dS/dpi^ {i = 1,2, 
and 4). They are conjugate to pj*. For coherent notation we also define 7^ = Q. We can 
easily find 

7:=y, (i = 1,2, and 4), (148) 



and the inverse is 



1' = ll 



7^ = Ai7i + A27* +Q + ^7* 



2m , (149) 
A 



4 4 
7 =7*- 

As for the momenta, from 

dS 

Pi = — = a,j-7J + Pi (150) 
(and the above relations for 7*) we find 

Pi = (an + ai3Ai)7^^ + {au + ai3A2)7^ + aisQ + (au + q;i3^)7* 
+ - AiP, 

P2 = ("12 + a23Al)7^^ + (a22 + 023^2)7* + "23(3 + ("24 + "23^)7* 
+ P2* - A2P, 

, A \ 1 / A N 2 ^ / 2m, 4 (151) 

P3 = ("13 + "33^1)7^ + (a23 + a33A2)7* + "33^ + ("34 + "33 — )7* 
+ P, 

Pa = (ai4 + a34Ai)7^^ + (024 + a34A2)7* + ^saQ + ("44 + 034^)7^ 
2m 

+ P4* yP- 

Note that the coefficients of 7* {i = 1, 2, and 4) in the above relations can be immediately 
explicitly expressed in terms of the background variables using the equations in (i) to (iii) . 
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The Gauge-Decoupled Form of Hamiltonian 

The Hamiltonian -ffeven for the new variables is obtained from 



Hp- 



-"even "r „ , 
Ot 



(152) 



where (d/dt) acts on the background variables only. The resulting Hamiltonian becomes 
of the form 

i^even = H^JeniQ, P) + JWl* + J2P2* + J4P4*, (153) 

where H^^^{Q,P) is the gauge-invariant Hamiltonian, and Jj = Ji{Q, P,jl,J^,jtjPi*) 
{i = 1,2, and 4) are the multiplier functions. 

Prom a straightforward computation we find the following gauge-invariant Hamilto- 
nian: 



^(llA^ + 4)7rig2 , (A^ + 2) 

H : 7ri7r2 



+ 



16 ■ 4 

A2(A2 - 1)% 5A4 



7r|g2 

5A2 



+ 



+ 1 Uo + 



A^(A2 - 1)91 



(A2 + 2)g2 



+ 4g2 A2 + 2 + "^^^^ 2 



q2 J (A2 + 2)g2 

"33 1 A^Q^ 



(154) 



033 2 ^ ^ ^VPQ 

2 2giyg^ 2yg^ 



where 



and 



$7 = g| + 2qiq2T, + 



A2 + 2 



2 ^^2 
^1 S , 



(155) 



2(A2-l)^2~^-^l + ^^2gr^(?2 



+ 



7A2 



+ 4 7r2 + 



5 
2 

2A2(A2-1)(;i 



7r2 92 



(A2 + 2) 



(156) 



+ (A2 + 2)g2 gi ■ 



As we noted earlier we have expressed aij in terms of 0:33 only, solving the equations 
cii = C2i = C4j = (z = 1 ~ 4) (except the "identities" C21 = C24 = C41 = 0). The function 
a33 is therefore now an arbitrary (prescribed) function of t. 
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As for the multiplier functions, we find 



NQ 



+ /„,^9 (A - 1)71-2 H + 77I"2g2 ' 



A 



J2 = - 



V2(A2 + 2) 
+ {'^2qiq2^ + 2^1033) E 



^712 + 2(A2 - I)7r2-^gig2 ^ 



A 



2 ATSgiP AV2(A2 + 2)N-fl 



A2 + 2 7r2v/^ (157) 



{X^^ 91 + 7r2g2) ^ + 



+ V2(A2 + 2)si, 



m iVQ 2mN^P 

Ja = — 7-(7r2 + 2g2a33)- 



+ 



A A 

2XmN-il 



7r2V92 



- (vTigi - 7r2g2) ^ H 

— msi + AS3. 

These functions provide the time derivatives of the gauge variables 7* {i = 1,2, and 4), 
i.e., 

ii = Ji\o (158) 
for each i, where |o stands for imposing the constraints pj* = 0. 
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The Wave Equation for the Even Perturbation 

Let us choose 033 so that it makes the Hamihonian a diagonal form. We therefore choose 



"33 = ~ 1)^2 ^ -^1 + ^^2^1 ^2 

Hl-^^j^-^ ^2.2 )WT2) 



(A2 + 2)g2 

to make Cgven = 0. Then, the Hamilton equations are given by 



^ = -11 16^^ + ^^-1-2 



A^(A^-l)^gi , 5A4 ^2 .\ 1 



4 y ' q2 y (A2 + 2)^2 

4^2 A^ + 2 J ^ 

It is easy (in principle) to find explicitly the wave equation for Q from the above 
Hamilton equations, but instead of writing down one that has very complicated form, let 
us present the wave equation for Q for the exact background solution ((T3|) . It is given by 



where 



X = 8m'^{t - 2A;)t+2 + 2 A^ m'^{2t - k) + 2A^(2t - 3A;) + ^-^^^ — 
Y^16mH^^-4k U (a^ ^ - s) + 2 A^ (162) 
Z = Am'^t+^ + 8m^ t- + 4 X"^ m'^ tt+'' + 2 X^ t_ t+ + A*^ 



8.2 Nongeneric (A = 0) perturbations 

In the following we work with a nongeneric even mode phase space So^rn and suppose that 
the perturbation canonical tensors are given by (See Ea. (|Hn|) ) 



lab = l^{El)ab+l^{E:i)ab, 



ab . . (163) 
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Note that dim^o.m = 4. 

As we have seen in Sec|lJ when A = we have two constraints DH^^^^^ ~ and 
T)Hi^~^^ ~ 0. (See Eas. H65|) and ((TOJ) Note, however, that T)H^^~^^ vanishes identically 
when m = 0, and thus we are left with a single nontrivial constraint I?i7{'^="'='') 0. We 
must therefore consider the two cases m = and m / separately. 



8.2.1 The A = and m ^ case: Pure gauge 

We have two (nontrivial) constraints 'DH(^x=o) ~ ^-ff«(A=o) ^ when A = and 
m ^ 0. Since the phase space So^m is four dimensional, these two constraints define 
4 — 2 = 2 dimensional constraint subspace and the gauge transformations these constraint 
functions generate on this subspace contract this subspace to a point: £o,m^o = {e}- In 
other words, the phase space <?o,m7^o is a space of pure gauge, and therefore we may not 
be interested in this space. 



8.2.2 The case A = and m = 0: Locally homogeneous perturbations 

In this case we have only one constraint I?i7{-^="^='') 0, so the dimension of the gauge- 
invariant phase space dimiSo,o is two(= 4—1 — 1). An explicit form of the constraint is 
obtained from Ea.(|H5|) by putting m = 0: 



2 2 

-p^(A=m=0) =_p^^^ _ (^^^ _ ^2 q^) _ j'^ 

+ 7 Ui + z ~ 0. 



(164) 



We do not have to perform computations from scratch to decouple gauge. As easily 
seen, our rule to reinterpret the computations presented in Sec l8.2l for the A = m = case 
is (1) first to put 7* = Pi = for i = 2, 4, (2) then take limit m ^ 0, (3) finally take limit 
A ^ 0. Keeping the order of steps (2) and (3) is necessary to drop correctly the terms like 
m/X. This order is particularly important when getting the correct wave equation, e.g., 
from En. fmTTl . 

We find that the gauge-invariant configuration variable becomes 

Q=-f^-^)y+^ (165) 

V^2 qij 

The wave equation can be obtained from the generic equation (|161|) by first putting m = 0, 
then taking the limit A — > 0, as noted. We have 

Q + ^Q = 0. (166) 

The solution is 

Q = Cit-^ + Ca, (167) 
where Ci and C2 are integration constants. 
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9 Gauge Properties of the Odd Perturbations 



In this section we discuss gauge issues for the odd perturbation system to obtain the 
perturbed spacetime metrics in a completed form. Remember that the system contains a 
gauge(-dependent) variable 7,, and the odd shift function s. Choosing a particular profile of 
the shift function s(t) is called a gauge choice for the odd perturbations. The evolution of 
the gauge variable 7=,, is determined by the evolution equation (the gauge equation) (jlUlf) . 
depending on a given gauge choice. The gauge equation also depends (parametrically) 
on the gauge-invariant variable Q, which is however determined independently of 7,, and 
s. By solving the gauge equation for a given shift function s we can express the gauge 
variable 7=,, in terms of Q, and thereby we can express the metric functions 7^ and 7^ in 
terms of Q. We see below how these functions are expressed in some gauge choices. 

We first need to obtain the map acting on the odd variables induced by diffeomorphisms 
for later use. 

Lemma 2 (Induced map) Let v{t) be an arbitrary function of time. Then, the one- 
parameter diffeomorphism generated by the odd vector 

Ya = v{t)Va (168) 

induces the following map /odd on the odd perturbation variables (5,7^,7^).' 



S + [ V V 

qi 



Iodd[v]: 75 ^ 75 + ^2(A2 + 2)v, 



(169) 



7^ ^ 7^ — mv, 



where s = s{t) is the odd shift function. 

Proof: This is a straightforward calculation of the induced map 

Sgab^Sgab + 2V^aYb). (170) 

Here 6gab is an odd perturbation spacetime metric, and Va is the covariant derivative 
associated with the background spacetime. A straightforward computation shows 

2V(^aYb) = 2(^v- ^^v^ (dt)(,H) + viV2{X^ + 2){E5)ab - m{EQ)ab), (171) 

from which we can immediately read off the induced map on 5Na, 7^, and 7^. ■ 
This induced map tells us that fixing the odd shift function s does not imply a complete 
fixing of the freedom of diffeomorphisms. We obtain the following. 

Proposition 3 (Pure gauge solution) Consider the odd perturbation equations \74\l 
with a given odd shift function s{t). Let Ci be a constant parameter. Then, the set of 
functions 

lU=Cm, ll^) = -C,^q,, (172) 

provides the one-parameter solution of the equations which is generated by spatial diffeo- 
morphisms. 
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Proof : The zero perturbation s = 7^ = 7^ = is a trivial solution of the perturbation 
equations. The diffeomorphism generated by the spatial vector H168|) induces the map 
(|169j) . which maps the zero solution to a pure gauge solution. We impose the condition 
that s = be retained, then this requires that we choose v(t) = Ciqi, where Ci is a 
constant. The pure gauge solution of the claim immediately follows after a redefinition of 
the constant Ci. This solution can be superposed with any other solution (5,7^,7^) and 
does not change the odd shift function s, from which the claim holds. ■ 

We can express the solution of the odd perturbation equations in terms of the gauge- 
invariant variable Q for a given gauge choice. While as we remarked, a gauge choice 
means a choice of the odd shift function s, it is often useful to express s with another 
gauge function. 

Proposition 4 (Gauge- variable-free general solution) Let Q be the gauge-invariant 
variable for the odd perturbation which evolves according to Eg. ilU4^ , and let 1) = ^(t) be 
an arbitrary function of time. Then, the solution of the odd perturbation equations i74\ ), 
up to diffeomorphisms, is given in terms of Q and by 

s = 0{Q,Q;^), 

75 = ^2(A2 + 2)i?Q, (173) 

76 = (1 _ ^^)Q^ 

where the function 0{Q,Q;'d) is defined by 

0.. ).(^-|.. 3^1) <,.(.- 3^^) (-) 

Proof: Remember that the gauge variable 7* is governed by the evolution equation 
'i* = J\o', see Eqs. (|100|) and (|lfllj) . For convenience we re-express P in terms of Q (and Q), 
using the Hamilton equations. Similarly, we express tti and tti (the background momentum 
variables) in terms of qi, q2 and their time-derivatives, using the background Hamilton 
equations. The gauge evolution equation 7^, = J\q then becomes 

'^*-{j^-fj^* = l ('^(^' ^) - ^) ' (1^^) 

where 0{Q, Q; 6) is the same as that defined in Ea. ()174() (with -d replaced by 9). Here, 6 is 
the mock shift function (the parameter function in defining the gauge variable 7=,,). From 
Eg. ipU)) . we can think of 7* as a function of 7^ (or similarly of 7^), 

7, = Q ^ ' (176) 

' ^ V2(A2 + 2)0 ^ ' 

Substituting this into the above equation we obtain the following evolution equation for 

^5 _ 91 5 ^ (Q\^ V2(A2 + 2)s. (177) 

91 u \qij 
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The homogeneous solution to this equation is given by "y^-^^^-^ = Ciqi, which is however 
found to be the pure gauge solution (see the previous proposition), so we drop it. To find 
the special solution we are interested in it is useful to transform like 7^ = f /qi- Using a 
new function we rewrite the shift function s as 

s = OiQ,Q;^) (178) 

so that the equation can explicitly be integrated. The solution for 7^ is then given by 
7^ = v^2(A2 + 2)i?Q. 76 follows from the definition of Q, 7*^ = Q - (1^/^/2)7^ ■ 

Note that there is no ^-dependence in Ea. H177|) . therefore 7^ and similarly 7^ do not 
depend on the mock shift 9. 

Remark The above solution can be viewed as the "zero-gauge" solution correspond- 
ing to the "gauge-fixing" 7=,, = 0. Remember, however, that we have freedom of choosing 
a function of time, the mock shift 6, to define the gauge variable 7=,,. Because of this 
arbitrariness the "gauge-fixing" by 7^, = does not have any invariant meaning. We can 
obtain the formula led by this observation. In fact, from Eq. (|175j ) the condition 7=,, = 
corresponds to taking the same relation H178|) (with 'd replaced by 9), and the same for- 
mula for 7^ and 7^ is easily found from Eq.(|9n|) and 7=,, = 0. Therefore the function is 
the same as the mock shift 9 under the condition 7=,, = 0. Nevertheless the above formula 
is valid regardless of this condition, and in such a general situation the two functions 
and 9 have nothing to do with each other; just parameterizes the odd shift s, while 9 
parameterizes the freedom of defining 7=,,. In particular, 1) can be zero, as opposed to 9. 

It is easy to express 7^ and 7^ directly in terms of the odd shift function s. To this we 
solve the relation s = 0{Q,Q;'d), which is a simple ODE for The solution is 

where we dropped the homogeneous solution ??(hom) = Ciqi/Q, which corresponds to the 
pure gauge solution. If we substitute this solution into Eqs. (|17:-{|) . we may obtain the 
explicit solution for a given odd shift function s. 

One of the particular gauge choices of our interest may be so-called the synchronous 
gauge, which is characterized by s = for the odd perturbations. 

Proposition 5 (Synchronous gauge) Let Q be the gauge-invariant variable for the odd 
perturbation which evolves according to Eg. \1U4\I - Then, the odd perturbed metric with 
vanishing shift is up to diffeomorphisms given by 



1 



V-2uq, (9. - [Q (Sl\dt 
\u J qi \ uy 

'^Q+.\J'^mdt. 

u J qi \ uJ 



(180) 
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Proof: This can be immediately obtained by setting s = in Eq. (|179) ): 



We perform an integration by parts then substitute the resulting -d = 'Sq into the general 
expression (|173() . ■ 

Formula H173|) provides the most general form of the odd metric functions expressed in 
terms of the gauge invariant variable Q. Let us confirm that they are certainly isometric 
to each other for any different choices of t?. 

Proposition 6 (Equivalence of the solutions) The odd perturbed spacetime metrics 
specified in Proposition^ with different i?(t) are all isometric to each other (up to the 
linear order). 

Proof: We show that the perturbed metric with an arbitrary in Proposition |1] is 
isometric to the one in Proposition which in turn implies that all the metrics are 
isometric to each other. Let 5^'^^ +e6g + 0{e^) be the perturbed spacetime metric that the 
perturbation metric in Proposition |3] defines. Let Ya = v{t)Va be an odd vector on this 
perturbed spacetime, and consider the diffeomorphism generated by this. From the map 
(fr69|l we find that 



7^ = V2(A2 + 2)i9Q^ v'2(A2 + 2)i}Q + y/2{X^ + 2)v 

= y^2{X^ + 2){^Q + v), ^^g2) 
7^ = Q — m-OQ Q — m'&Q — mv 
= Q-m('dQ + v), 

from which we see that the transformation by the diffeomorphism is equivalent to trans- 
forming 

= ^Q + v. (183) 

On the other hand to make 6Na = 0, it is found from the map (|169() that we need to 
choose v{t) to satisfy 

v-^v + 0{Q,Q;^) = 0. (184) 
We can easily solve this first order ODE, and obtain 

„(0 = -,J^(«lMd, (185) 

J qi 

With this V, we have 

^> =^Q + y = ^Q_q^ [ 9iR2Ql3.dt. (186) 

J qi 

What we need to confirm is that this 'd' does not depend upon the original choice of and 
coincides with 'do defined in Ea. (|181() . To this we first differentiate the above equation 
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(after an appropriate rearrangement) to delete the integral. Then, we have the following 
ODE for 1?': 

(187) 




Substituting the definition (|174() of 0{Q,Q;i9), we obtain 




^' + f ^ - - I ^' = Tf^- ( ^ - - 1 • (188) 



Q qi A2 + 2 n 



Now, we can see that the apparent dependence upon the original choice of has disap- 
peared. Moreover this equation is the same as 0{Q,Q;'&') = 0, which implies i}' = t?o- 
Hence, all the perturbed metrics with different i) are isometric to the synchronous per- 
turbed metric given in Corollary |S1 ■ 
Finally, we introduce a particularly useful gauge choice in the present perturbation 
system. With this choice the perturbed metric does not contain, in contrast to the syn- 
chronous form given in Corollary |SJ the time-derivative of Q and can be written in the 
simple multiplicative form 6g = Q x G, where G is a symmetric tensor that depends only 
on the background variables. 

Proposition 7 (Standard gauge) Let Q be the gauge-invariant variable for the odd 
perturbation which evolves according to Eg. ilU4^ . When the perturbation of the spacetime 
metric is written simply in the product form 6gab = QG^^j^^^ for a symmetric tensor G^^j^'^^ 
which is independent of Q, the tensor G^^'^*^^ is given by 

= V2.'-iE,)ab + 'AE,)ab + 2-^ r-^)Viadt,y (189) 
u u y + 1 \ u J ^ ' 

Proof: Let us choose the function •& in formula (|17.3j) so as to eliminate the Q- 
dependence of the perturbed metric. When m 7^ 0, we can eliminate the Q-dependence in 
0{Q,Q]^) by choosing (see Eq. (|T7H) ) 



and in this case we have 



u 



5 = V2v^Q, (191) 
u 



7 

u 

This completes the proof for the case m 7^ 0. Since it is possible to take limit m — > 
in the above expression, the same formula is correct in the case m = 0, as well. [In fact, 
when m = the solution of (| 177(1 is given by 



1' 



V'2(A2 + 2)gi / -dt (192) 
J Qi 
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for a given s. From the assumption, we must have s = fQ for a function / which is 
independent of Q, but in this case the above 7^ is not of the assumed form unless / = 0, 
which reproduces Eqs. (|191|) for m = 0.] ■ 

This gauge choice is similar to so-called the longitudinal gauge known in the pertur- 
bation theory of isotropic and homogeneous spacetimes. We call this gauge choice the 
standard gauge for the odd perturbation, due to its simplicity. 

Other possible choices of gauge which are of some interest would be those by 7^ = 
("solid base" gauge), and 7^ = ("orthogonality-preserving" or simply "orthogonal" 
gauge), which, respectively, correspond to = and -& = l/m. ("Orthogonality" here 
refers to the one between the base T,g and fiber S^.) 

10 Gauge Properties of the Even Perturbations 

Let us first see how diffeomorphisms act on the even perturbations. 

Lemma 8 (Induced map) Let ai (z = 0, l,3j he arbitrary functions of time. Then, the 
one-parameter difjeomorphism generated by an even vector 

Ya = aoS{dt)a + aiSa + a3Sa^a (193) 
induces the following map /even on the even perturbation variables: 

si ^ si + di — q^^qiai + Xoq, 
•S3 ^ S3 + «3 - q2^<i2ai - mao, 

r r 1 1 1 ^1 , (194) 

Uwen[ao,ai,a3\ : 7 ^7 -j^aQ-Xai, ^ ' 



' -2 + ^2(A2 + 2)ai, 



7^7 



3 3 '72 
7^7- + 2ma3, 

7'* ^ 7"^ - ma\ + Aas, 

where sq, s\, and S3 are, respectively, the even lapse function, even base shift function, 
and even fiber shift function. 

Proof: This is a straightforward computation of 2V(aYi,). For our background metric 
gab = -N'^{dt)a{dt)b + qab, it holds that 

/ 1 . ^'''^ 



-{dt)a{dt)b 2^2- 
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Also, for a spatial vector Aa such that Aai^dt)"" = 0, we find 

VaAb = DaAb - q"^Acqdia{dt)b). (196) 
Using these relations we can find 

2V(„n) = 2S ^do + ^ao^ {dt)a{dt)b 

+ 2 (di - q^^qiai + Aoq) (dt)(a5'fe) 

+ 2 {as- q2^q2ai - mao) S (dt) (^^(^^ b) (197) 
- {qi{Ei)ab + q2{Ez)ab) 



ai iX{Ei)ab - ^2(A2 + 2)(^2)a6 + m{Ei] 



ab 

+ as {2m{Es)ab + KEi)ab) , 

from which we can read off the map claimed. ■ 
It is straightforward to confirm that the even gauge-invariant perturbation variable 
Q = — Ai7^ — A37^ + 7^ — ^7^ is certainly invariant under this map, Q ^ Q. (To 
confirm this we may need the relation Ai = q2qi^ , which is obtained from the definition 
(|144|) using the background Einstein equations <\l2\i .) 

As in the odd case let us find the pure gauge solution that can be superposed on any 
solution with given even lapse and shift functions. 

Proposition 9 (Pure gauge solution) Consider the even perturbation equations with 
a given even lapse function sq and even shift functions si and S3. Let Cq, Ci, and C2 be 
constants and let 

Co 

"° = ^' 

/di 
iV^ + Ciqi, (198) 

as = mq2Co [ + C2q2- 



Then, the set of functions 



Nq2 



1 ?i X 



^(d) = -^ao + 2ma3, 
T^d) = -max + Aas, 

represents the three parameter solution of the equations which is generated by diffeomor- 
phisms. 
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Proof: Use the same logic of Proposition |21 Consider the diffeomorphisms generated 
by the even vector of the form H193|) . From Eas. H194|) we find that to keep 5Na = 

ai - qi^qiai = -Aao, 



_i . (200) 
03 - 92 9201 = mao 



should hold. The solution is easily found, provided that oq is given: 

ai = -Xqi / h Cigi, 



; % (201) 
03 = mq2 / h 6292- 

On the other hand, we also find that to keep 5N = 

N 

do + ^ao = 0. (202) 

The solution is 

ao = CqN-^. (203) 

The claim follows by substituting this into the above ai and 03, and applying the corre- 
sponding /even map to the trivial (i.e., zero) solution. ■ 

As in the odd case, we can find the longitudinal gauge-like solution, where the gauge- 
invariant variable Q appears in a simple multiplicative form. 

Proposition 10 (Standard gauge) Let Q be the gauge-invariant (configuration) vari- 
able for the even perturbation which evolves according to Eqs. Ufif^) . When the perturbation 
of the spacetime metric is written in the product form 6gab = QG^ab^^^ ^ symmetric 
tensor G^'j^™"^ which is independent of Q and depend only on background quantities, the 
tensor G^^j^^^^ is given by 

~,(even) 1 



ab 



qiq2{Ei)ab + 2qtA2{E2)ab 

TTl 

- 92(giAi + q2)iE3)ab + 2ygig2(^4)ab 

+ A2T2 " + 392) S{dt)a{dt)t 

o/V^ / \2 \ 

+ f Ti + 3^2^T2 j Sidt)aidt)i, 

- A2T2 + T3) (AS(„(dt),) + mSa\,idt),)) 
Am 



(204) 



where 



-^'^iSa (a{dt)b), 
Ti-m. T.-f4?Y. T3^4P. (205) 



and symbols Ai, A2, S, and Q are those defined in Eqs. \14'^ , j j^5| ), and il55\) . 
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Sketch of the derivation of i204\ )- The assumption of the gauge leads us to put 7* = 
(i = 1,2,4), where i?j are functions of the background quantities. (7^ can be determined 
from the definition of Q.) From Eqs. (|148j) . this also means 7* = -i^jQ {i = 1,2,4). We 
substitute these relations into the gauge equations H158() . We then eliminate P and rewrite 
the gauge equations in terms of Q and Q by using (one of) the Hamilton equations (|160j) . 
and determine the functions 'di so that the coefficients of Q in those equations vanish. 
Those equations are simple algebraic equations which can be immediately solved, and the 
solution is written in terms of the background quantities only. We have 



^1 



^4 



(M2 

' n ' 

2gfA2 
2m gig2 



(206) 



From these 
functions s,- 



7* {i = 1,2,4), and therefore 7^ are determined. The even lapse and shift 
(z = 0, 1, 3) are determined from the gauge equations with the above choice of 
i^i. They are also algebraic (linear) equations for Si, so it is easy to solve, and the solution 
automatically satisfies the gauge assumption. We have 



So 



si 



S3 




(207) 



(We have used the background constraint W = to simplify the result, e.g., eliminating 
one of the background variables vri in sq.) From these we obtain (|2U4() . ■ 
Due to the complexity (or lengthiness) of the even gauge equations it may be im- 
practical to write down a general solution for an arbitrary gauge. Instead, we can derive 
any solution with a desired gauge by acting the nicip leven tli6 a,bovG stciTiddTd gauge 
solution. In this case, the functions which are considered free are the three functions 
(i = 0,l,3). 



II Summary and Discussions 

A brief summary first. We have seen that mode functions on a closed orthogonal Bianchi 

III space are specified by two kinds of eigenvalues; the fiber eigenvalue m and the base 
eigenvalue A. Using the metrics and covariant derivatives on the base and the fiber 
S^, we have constructed the "even" mode vectors and tensors from the mode functions. 
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Using the area two-form Sab on the base in addition to the same metrics and covariant 
derivatives, we have constructed the "odd" mode vectors and tensors. The "harmonic" 
mode vectors and tensors are associated with the (Hodge) harmonic vectors on the base 
and the mode functions on the fiber. The "tt" mode tensors are associated with the tt 
tensors on the base and the mode functions on the fiber. This way the mode quantities 
can naturally be categorized into the even, odd, harmonic, and tt ones. Perturbations 
of vacuum solutions associated with different kinds (i.e., even, odd, • • • ) and/or different 
eigenmodes are all independent of each other, in other words, they all decouple, as stated 
in Theorem^ This makes it possible to treat them separately. 

We have defined mode tensors {Ei)ab (i = 1 ~ 9) for a given set of m and A (i = 1 ~ 6) 
or for a given m (i = 7 ~ 9), which are eigentensors of the Laplacian and are L^-orthogonal 
to each other with respect to the standard L^-norm. To separate variables we expand the 
first variation of the spatial metric in terms of the mode tensors, resulting in getting 
the coefficients 7* = 7*(t) as the fundamental perturbation configuration variables. The 
canonical conjugates can also be expanded by the same mode tensors and we obtain 
Pi = Pii't) ^ the conjugates of 7*. The Hamiltonians for these perturbation variables are 
provided in Eg. ((7^ (odd part) and H125() (even). That the harmonic part can be obtained 
from the odd one by formally taking limit A ^ was pointed out in SeclHl The Hamilton 
equations for the tt part are given in Eqs. (|114j) . 

A set of Hamilton (i.e., evolution) equations for perturbation variables is generally 
a system of complicatedly coupled linear ODEs which depend on choice of gauge, i.e., 
depend on choice of variations of lapse function 6N and shift vector 5Na- One of our 
main tasks in our course of analysis is to decouple the system into gauge-dependent and 
independent parts for each one of the even and odd perturbation systems. (The harmonic 
one can be thought of as a limit of the odd system, as already remarked. Since the tt 
system is not constrained, the fundamental variables for the tt system are themselves 
gauge-invariant.) One of the advantages of employing Hamiltonian formalism is that this 
task amounts to finding a certain kind of canonical transformation. 

We have used the method of generating function to find such desired canonical trans- 
formations. This systematic method was crucial for our achievement of the task. The 
decoupled-form Hamiltonians are given in Eq. with and UlUUf) (odd) and Eq. H153|) 
with (|154|) and (|157j) (even) . The wave equations for gauge- invariant variables Q (the defi- 
nition of Q is different depending on the type of perturbations), or the Hamilton equations 
for Q and its conjugate P are given in Eq. (|lfl4j) (odd) and (|16fl)) (even). These evolution 
equations for the gauge-invariant variables are the "heart" of the perturbation dynamics. 

To recover an entire perturbed spacetime metric, however, it is unavoidable to choose 
a particular gauge. As a useful choice, we have defined the standard gauge, in which the 
gauge-invariant variable Q appears in the perturbed metric in the simplest manner where 
the metric is expressed in terms of Q only and free from Q or P. The forms of metrics in 
this gauge are given in Propositions 171 and ITUl 

We have thus resolved all the fundamental issues concerning the technical aspects of the 
perturbations of closed orthogonal Bianchi HI vacuum spacetimes, including separation 
of variables, decoupling into gauge-dependent and independent parts, derivation of wave 
equations, and choice of gauge. 
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Our primary interest concerning properties of the perturbations of the Bianchi III 
spacetime may be their future asymptotic behaviors. For this not only do we need to 
perform a further careful analysis of the wave equations but we need to clarify how to 
subtract the effect of the anisotropic expansion of the background and thereby get the 
"net" effect of the perturbations. This nontrivial task is indispensable to see under what 
condition the system should be called "stable," or to say at what rate the perturbations 
decay or grow. These subjects will be discussed in a separate paper. In the remaining 
part, let us discuss the scalar and electromagnetic field equations and the tt perturbation 
equations, which are simpler than the other perturbation equations. 

In Appendix El we have dealt with the wave equation for massless scalar fields. When 
the background is Minkowski (with a slicing that is Bianchi III), we can obtain the exact 
solutions for the (mode-decomposed) wave equations for all possible eigenvalues. The 
solutions tell us that the fields with m 7^ all decay as t ^ cxd at the same rate (See 
Eq. (|250j) ). while the C/(l)-symmetric (i.e., m = 0) fields decay at several different rates 
depending on the value of A (See Eq.^E^). It is however true that all modes except 
the zero mode are decaying. It seems therefore plausible that all the modes are decaying 
(except the zero mode) for any general orthogonal Bianchi III backgrounds, too. 

The (source-free) electromagnetic field equation has been dealt with in the subse- 
quent appendix. Because of the vector- like nature of the electromagnetic fields, the mode- 
decomposed wave equations belong to either the even type or odd type (or the harmonic 
type, which can be obtained as the limit A — > of the odd type). Interestingly enough, 
the wave equations for both even and odd fields coincide with those for the scalar field 
(after a normalization if necessary) when the background is Minkowski. They are therefore 
decaying again on this background (if we think of the normalizations involved natural) . 

Among the gravitational perturbations the simplest are the tt perturbations. This 
kind of perturbation is, as noted in SecjTl the only exceptional case among the four kinds 
where it is clear how the variable should be normalized to see the stability. We have 
seen in the same section that the wave equations for the (normalized) tt perturbation 
variable, and therefore their asymptotic behaviors, are exactly the same as those for the 
scalar fields (for A = 0). The background is therefore expected to be stable against the 
tt perturbations. 

As noted, detailed studies of dynamical properties of the perturbations will be reported 
in a separate paper pi|. 

A Spatially Closed Bianchi III Spacetimes 

In this section we discuss how we can construct the spatially closed Bianchi III spacetimes. 
The so-called Teichmiiller parameters are the key ingredient to understand the compact- 
ifications involved. The Teichmiiller space of a manifold M modeled on a homogeneous 
space X is the quotient space J^/Vq, where M. = M.{M,X) is the space of the locally 
homogeneous metrics on M which are locally isometric to X, and T>q = T>q{M) stands for 
the diffeomorphisms of M that are connected to the identity. Teichmiiller parameters are 
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defined as coordinates in a Teichmiiller space and they are interpreted as the config- 
uration variables that emerge as a result of compactification. We note that Teichmiiller 
parameters can also be understood as the parameters in the covering group acting on X. In 
the following we first show how we can compactify Bianchi III homogeneous manifolds by 
explicitly finding covering groups acting on them, and then generalize the compactification 
to ones for spacetimes. 

Let = {Y{?,hab) be the standard hyperbolic plane, and = (R, /q;,) be the stan- 
dard Euclidean line. We introduce the standard coordinate systems with which these 
metrics are expressed^'^ 

hi^dx'dx^ = '^^AfiL (208) 

and 

kjdx'dx^ = dz^. (209) 

The coordinates are defined in the region x £ (— oo,cxd), y £ (0,cxd), and z G (— oo,cxd), 
and they span the whole and E^. 

As well known can be compactified to a higher genus surface {g > 2) by the left- 
action of a discrete subgroup F''^ of the orientation-preserving isometry group Isom+i?^: 

tth: H^^ (S3, hab) = i/Vr^- (210) 

The topology of the closed surface can be specified by its fundamental group 7ri($]g), 
which is an infinite group generated by 2g generators Oi and /3j (i = 1 ~ g) which satisfy 
a single relation [ai,/3i] • • • [ag,f3g] = 1. That is, in the standard notation, 

'^li^g) = {ai,--- ,ag,Pi,--- , Pg\[ai, Pi] ■ ■ ■ [ag, Pg] = 1), (211) 

where [ai,/3i] = aiPta'^^fS^^ is the commutator of group elements. The covering group 
is a representation of this infinite group into Isom+Z/'^. 

The group Isom+i/^ is isomorphic to the projective real special linear group of second 
rank, Isom^H^ ~ PSL(2, R). Suppose we embed the generators Oi and /3j into PSL(2,R) 
so that they satisfy the relation for TTi(T,g). (That is, we identify Oi and /3j with elements 
in PSL(2, R) so that they satisfy the relation for 7ri(Sg) with respect to the multiplications 
in PSL(2,R).) However, we do not need an explicit representation of such and (5i G 
PSL(2,R). We just note they generate F''^, which we denote 

F^ = {ai,--- (212) 

where a^, /Si e PSL(2,R) (i = 1 ~ g). 

^^Beware that several different definitions of "Teichmiiller space" are used in the literature. In general, 
they are equivalent for two dimensional surfaces, but not for higher dimensional manifolds, including 
three-dimensional ones. Ours follow [113 11^113 HOI HI]- 

^''The hyperbolic plane model with this metric defined in the region y > is known as the upper-half 
plane model, and the plane is often denoted like R+ to emphasize it is defined on the upper-half plane. In 
this section however we understand that is homeomorphic to and do not distinguish those two for 
notational simplicity. 
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Note that, since PSL(2, R) is three dimensional, we need (at most) 2g x 3 = 6g 
parameters to represent . However, since the generators must satisfy the relation 
[ai,/9i] • • • [ag,f3g] = 1, the number of independent parameters reduces by three, because 
this relation has three components. Furthermore, remember the fact that different 
and r^' can result in the same closed surface, i.e., /T^ and H^/F^' are isometric to 
each other, if = a^^ o T^' o a for an a G IsouiH^. (This map is called conjugation. 
See, e.g., [IJ.) Since Isom/7^ is three dimensional, these conjugation maps reduce the 
number at most by another three, and we know they do reduce it by the maximal number, 
hence resulting in 6(7 — 6 independent parameters in . This count corresponds to the 
well-known fact (e.g.p^j) that a higher genus surface with genus g > 2 has 65 — 6 
Teichmiiller parameters. 

Let us generalize the covering group of the closed surface to that of the trivial 
S'^-fiber bundle M ~ x S"^. The topology of this closed manifold can be uniquely 
specified by its fundamental group, which can be represented by 

7ri(M) = (ai,- • • ,ag,/3i, - • • , (3g,tp\[ai, (ii] - ■ ■ [ag,(5g\ = = [(3i,ip\ = 1) , (213) 

where the new generator tp corresponds to the loop along the S'^-fiber. (This expression 
of TTi follows from the general theory of S'^-bundles over a surface or more general Seifert 
fiber spaces. See, e.g., |H1. ) 

We consider, as the universal cover, the direct product x E'^ with the rescaled 
metric 

Qab = qiKb + q2iab, (2i4) 

where qi and q2 are constants. The isometry group of this manifold is the direct product 
of those of and E^, i.e., Isomff^ x i?^ ~ Isomiif-^ x Isomi?^, where the orientation- 
preserving component of Isomi?^, Isom+i?^, is the one-dimensional translation group. So, 
an element a G Isomo//^ x i?^ ~ lsom.-^-H^ x Isom_|_£'^ can be^^ expressed as 

(215) 

using a G Isom+Z/^ and s G Isom+i?"^. This acts on points x = {x,y,z) G M, viewing 
x = (x,y) X {z} G X E^. The multiplication for two isometrics a, a' G Isomg-ff^ x E^ 
is defined by 

where aa' is understood as the multiplication in PSL(2,R). 
Proposition 11 Let 

r^ = {ai,--- ,ag,/3i,--- (217) 

^''The group Isomo//^ x ~ Isom+H^ x Isom+iJ^ is the identity component of Isom+i/^ x , the 
orientation- preserving isometry group of x . We focus on this component to reahze a compactification. 
To generahze the argument to the full Isom+iif^ x E^ , one may need to consider the orientation-reversing 
isometrics of and E^ , as well. 
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be the covering group for {Tig, hah) acting on . Let M be the product T,g x and let 
{M,qab) be a locally homogeneous manifold modeled on x = (Rp,qab), where ijab is 
the rescaled metric \214\). Let us represent {M,qab) with the following covering map 



^ : (R3, qab) ^ (M, qab) = (R^ gab)/r. (218) 
Then the covering group T C Isomo-ff^ x can be represented as 



V, h'"\vg j'l u;i y''"'UJ'U 



(219) 



where G R, u^j G R and (/? > are arbitrary constants. Here, the 8g — 5 independent 
parameters in T, i.e, 6g — 6 parameters in plus the 2g + 1 parameters Vi, Wi and ip, 
span the Teichmiiller space of M modeled on {Il^,qab)- 

Proof: Let us denote the generators of F as in 

T = {ai,--- ,aaJir-- Jg,^}, (220) 

where di,Pi,ip G Isomo-ff^ x E^. From the relation [ai,/3i] • • • [ag,f3g] = 1, it is obvious 
that the Isom+i/^-part of Si and /3j should satisfy the unarrowed same relation, from 
which we have 

"•=(::)' 

where Vi and Wi are real constants, and ai and f3i are the generators of F-'^. Conversely, if 
Si and Pi are given in this form, from the fact that Isom-|_£'^ is abelian, we see 

[aiJi]---[SgJg]= (^l^ = 1, (222) 

that is, the 2g generators (|22H) satisfy the required relation for tti(M) for arbitrary Vi and 

Wi. 

For an isometry a G Isomo-ff^ x E^ of the form 



(223) 



the conjugation map acts on Si (and similarly on /5j) as 



a^oSioa=l ^ j . (224) 

That is, the conjugation is effective only on the Isom+iJ^ part. From the discussion on 
F^, this fact merely implies that the 2g generators Si and jSi contain 65 — 6 independent 
parameters. 
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It is also easy to check that to meet the remaining commutativity relations [cSj , ip\ = 
\j3i, Lp\ = 1 the S^-fiber generator should be of the form 




(225) 



where > is an arbitrary positive parameter. (99 cannot be zero for M not to degenerate. 
Moreover, (p = (1, ip) is conjugate to (p- = (1, —ip) by the reflection isometry r : (x, y, z) — > 
(x,?/, —z). So, it is enough to consider positive ip.) The conjugation maps with respect to 
the isometries connected to the identity do not affect this generator. This completes the 
proof. ■ 
The Teichmiiller space of M was first discussed in jlUl I16j . 

Figure n pictures a fundamental domain of a compactified Bianchi III manifold and its 
projection onto the ff^-plane of z = 0. 




y 



Figure 1: A fundamental domain of a compactified Bianchi III manifold (the decahedron 
in the figure) and its projection onto a i/^-plane (the octagon in the x-y plane). (For 
convenience of presentation, the x-y coordinates in this figure are assumed to be those of 
the Poincare model rather than the upper half plane model.) 

Note that the universal cover x E"^ has two natural structures; the foliation by E"^^ 
and the foliation by . It is easy to see that the foliation by -E^, which is equivalent 
to the foliation by our z-axes, descends to a foliation by after the compactification 
we discussed (cf. Theorem 4.13 in |[12i)- In contrast, the above proposition tells us that 
the if ^-foliation does not descend to a S^-foliation unless vi = Wi = Q for all i. (Note 
that the iJ'^-foliation can be expressed in terms of the coordinates as the foliation by 
z = constant surfaces. So, the condition that the ii^-foliation descends to a S^-foliation 
is that each z = constant surface is compactified to S^. This does not happen unless 
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Vi = Wi = 0, as seen from the previous proposition.) We call these 2g parameters Vi and 
Wi the distortion parameters. It may be natural to consider the subclass where the 
inheritance //^-foliation Sg-foliation holds: 

Definition 12 When the distortion parameters all vanish, i.e., when the covering group 
r in Proposition \11\ is of the form 

{(oO--(oO'(o)-'(oO'(:)}' 

we say that the resulting closed manifold (M, Qab) is orthogonal. 

Note that the closed hyperbolic surface does not exist (or, is not embedded) 

in {M,qab), unless {M,qab) is orthogonal. 

When (M, qab) is orthogonal, the manifold admits both the foliation by (Sg, qihab) and 
foliation (or fibration) by {S^ ,q2lab)- The metric qab which is induced from the universal 
cover one qab can be written in the form 

Qab = qiKb + q2lab- (227) 

A tangent vector = v^dx + v'^dy + v^dz on p G M can be uniquely decomposed as 

v'' = v% + v''s, (228) 

where v'^ = v^dx + v^dy is the part tangent to the (Sg, gi/iab)-leaf to which the point p 
belongs, and v'g = v^dz is the part tangent to the {S^, g2^a6)-fiber to which p belongs, vfj 
and v'g are orthogonal to each other with respect to the metric qab- We can also check 
from the above form of the metric that the covariant derivative operator Da associated 
with qab can be written in the form 

Da = ba + cj\xz, (229) 

where Da is the covariant derivative operator associated with hab and X3 = dz is the 
derivative operator along the circle fiber. It is easy from this formula to see that the 
Laplacian = q"-^DaDb with respect to qab is given by the sum 

^, = q^^ ^h+q2\x^)\ (230) 

where A/^ = h°'^DaDb is the Laplacian with respect to hab- Therefore the product of 
eigenfunctions for {T,g,hab) and for {S^,lab) becomes an eigenfunction for {M,qab). 

Our next step is to generalize our compactification to one for spacetimes. Let Y = 
(R^,gaf,) be the universal covering spacetime of a spatially closed Bianchi III spacetime 

^^We do not use the word "twist" as in 9 to avoid confusions with those in the so-called pants- 
decompositions of hyperbolic higher genus surfaces. 

^''It is apparent that even if {M,qab) is not orthogonal there exist a foliation of M by Eg, since M is 
(topologically) the product Eg x S^. The induced metric /i^j, on each leaf of such a foliation is however 
not hyperbolic, h'^i, ^ qihab {qi is a constant). 



53 



Y = (M X H^Qgij). Let t be a time function of (R^,gafe) defined in such a way that each 
t = constant surface coincides with (R^^qab), where the scale factors qi and q2 in ijab are 
regarded as functions of t. Then, we assume that the universal covering spacetime metric 
is of the synchronous form 

gab = -N^{dt)a{dt)b + qab 

= -N {dt)a{dt)b + qihab + q2lab, 

where = N{t) and {dt)a = ^ at- Apparently, IsomY = IsomiJ^ x Isomi?^. Isomy is also 
equivalent to the isometry group of the spatial manifold, Isoml" = IsomM = Isom(M, qq^)- 
In terminology adopted in j2Ull21j . this fact can be rephrased as EsomM = IsomM, where 
EsomM is the extendible isometry group of (R^,gab)- In general, n = dim IsomM — 
dim EsomM gives the degrees of freedom of choosing the velocities of the Teichmiiller 
parameters of M at an initial surface t = constant |2fll I21j . In the present case where 
n = however there is no such freedom, and the compactification of the spatial manifold 
straightforwardly generalizes to the compactification of the spacetime manifold, i.e., our 
spatially closed spacetime solution can be represented as 

(MxR,g,,) = (R^g,fc)/^, (232) 

where F is of the form H219|) and acts on R^ as 

r • (t,x) = (t,rx). (233) 

gab may be assumed to satisfy Einstein's equation on the universal cover R"^. We can write 
the spacetime metric gab as 

gab = -N^{dt)a{dt)b + qab (234) 

in both cases where (M, qab) is orthogonal and non-orthogonal. It is however only for the 
orthogonal case that the spatial metric qab can be decomposed as Eq. (|227|) . 

The vacuum universal cover metrics gab form a one parameter family up to isometry. |17| 
On the other hand, T has 8(7 — 5 parameters, which are related to the Teichmiiller pa- 
rameters of the spatial section. The general solution is therefore a 8(7 — 4(= 8(7 — 5 + 1) 
parameter solution. The orthogonal solution, where T = Forth) has only 6(7 — 4 param- 
eters. However, even if the background solution is of orthogonal type, i.e., even if the 
distortion parameters of the background solution are all zero, it is worth remarking that 
the distortion parameters can be perturbed off the zeros (through the zero mode harmonic 
perturbations). 

^'^The isometry group of a spacetime with time function t has the subgroup whose action retains the 
foHation hy t = constant. This subgroup can also be identified with a subgroup of the spatial manifold 
specified by t = constant, and called the extendible isometry group of this spatial manifold 1201 1211 . Note 
that in general to compactify a spatially homogeneous spacetime y to a spatially locally homogeneous 
spacetime Y with spatial manifold M, we need to embed iii(M) not into IsomAf but into EsomAf. (This 
is apparent if we think of EsomAf C Isomy.) 
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B Wave Equation for the Scalar Field 

Consider the scalar wave equation 

g^^VaVfc^- = 0. (235) 

This equation can be derived from the action / = (—1/2) J g'^''(Va^')(V5^')(i//g, from 
which we can read the following Lagrangian for our background 

L = -I^I^Vq^ Jj-^(^)' + q''\Da^){D,^))di,,,. (236) 
The momentum p\ji conjugate to ^ is 



d^i N 



so the Hamiltonian is given by 

H = pq,4/ - L 



= ^ / + qi^2q''\Da^){D,^)) dfi,,. (238) 

Let us expand the wave function using the mode function S = Sa,™ as 

^ = ^PS, p^=p^S, (239) 

where ■0 = ^{t) and p^ = Ptp{t) are functions of time, since each mode of this form 
decouples from the others. Substituting this form into the above Hamiltonian and per- 
forming integration by parts (and dropping the constant / S'^dfig^ = 1), we have a reduced 
Hamiltonian 

H = ^ + ,i^(A2gr^ + m'q,')^pA . (240) 

The Hamilton equations derived from this are 

(241) 

or eliminating p^ we have a wave equation for ip, 

+ + NHX'ql' + m'q2')ij = 0. (242) 

yqi ^q2 y 

Upon substituting the exact solution (|15|) we obtain an explicit equation 

+ f A^^ + m^^j ^^ = 0. (243) 
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Another practical form is obtained by transforming {t-\.t- /tY^'^if): 



C + \C + Zsu,C = 0, (244) 



where the coefficient Zg^ of C, is 



(245) 

Solutions for the Minkowski background (A; = 0) 

Let us consider the case when k = (the Minkowski background). In this case the 
wave equation ()244() becomes Bessel's equation 

C + ^C+(m^-^)c = 0. (246) 

The original field and ( are related through tp = C/Vi- We can express the general 
solution to the above equation for ( with Bessel functions (when m ^ 0): 

C = Ci Ji^i34A^(mt) + C2iVi^/x34A^(mt), (m / 0) (247) 

where Ci and C2 are integration constants. J and N are, respectively, the Bessel function 
of the first kind and that of the second kind. It is worth noting that when A = this 
solution can be written only with elementary functions: 



( = CiW^-sinmt- C2A/^-cosmt. (m / 0, A = 0) (248) 

The asymptotic form of the solution (|247|) as t ^ cxd is of our particular interest, which 
is, as found in any appropriate formula table, given by 



. ^ / 2 / (1 + Vl-4A2)^\ ^ 12. / (1 + Vl -4A2)7r 
C~ Cn/ cos mt-^ — +C2W sm mt - ^ — ^— 



TTUit \ 4 / V irmt 

^ (249) 

Hence, in the range A^ < 1/4 the changes of A result in changes of phase in the asymptotic 
form, while for A^ > 1/4 they change the amplitude for a fixed t (since, e.g., 5R(cos(a+6«)) = 
cosh 6 cos a). The decaying rate (of the amplitude) of the field = C/V^ is, however, 
universal, which is given by 

|V'|=0(^). (t^oo,m/ 0) (250) 

The equation for m = should be considered separately, but finding the exact solution 
is not difficult, which is given by 

+ C2t-5v^i^, (m = 0, A2 < 1/4) 

Ci + C2 log t, (m = 0, A2 = 1/4) (251) 

Ci cos fW4A2 - llogt) +C72sinfiV4A2-llogt) , (m = 0, A2 > 1/4) 
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where Ci and C2 are integration constants. The decaying rate is therefore 

'0(t5(Vi-4A2-i))^ (m = 0,A2<l/4) 

O(^), (rri = 0,A2 = l/4) (252) 

O(^). (m = 0,A2>l/4) 

A critical point appears at A'^ = 1/4 in the present m = case, as opposed to the m 7^ 
generic cases. (We remark that the critical eigenstate for A^ = 1/4 does not necessarily 
exit. Also A such that A^ < 1/4 (so-called "small eigenvalue") may not exist. The spectrum 
varies depending upon the compactification of the background.) 

C Electromagnetic equations 

Consider the source-free Maxwell equation 

g'^'VeF.b = 0. (253) 
Here, the electromagnetic tensor Fab is obtained from the vector potential Aa through 

Fab = daAb - dbAa. (254) 
In local coordinates Eq. (|253j) reduces to 

= -La,(^g^"g'^^F„^) 



^'[^''f'\ '''g''F,o + ao(g^^g°°F,o) + -a.(/..g'^"g^^F«,). (255) 



The odd part 

For the odd part we can put 

Aa = a{t)Va, (256) 

where Va = CmVa = Cm^a^DbS. Equation (|255() is found to be trivial for fi = 0, and ^ = 3, 
so we are left with /u = 1, 2. In this case, Eq. (|255|) reduces to 

, _ao(iV.i V^)g^^gOO^^^ + 5o(g^^g°°F^o) 



^'^i^ (257) 

+ — 5z5(/i^g^^g^^Fsc) + 53(g^^g=^3^B3). 



^^Greek indices run to 4, and small Latin ones 1 to 3. g is the determinant of g^^. 
^^Capital Latin indices run 1 and 2. Any hatted object is lowered and raised by the standard hyperbolic 
metric Hab and its inverse h^^ . 
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Using 



and 



Fbo = -aCmVB, (258) 
Fb3 = maCrnVs, (259) 
Fbc = acmiDBVc - DcVb), (260) 



A/, Va = -(A^ + 1)Va, (261) 
b'^DAVB = -Va, (262) 



we obtain the following wave equation: 



" + ( - S I « + + m%^)a = 0. (263) 



2q2 

Upon substituting the exact solution we have 



2/c . A 9 1 9 ^- 

t+t_ V t+t- 



d + a+ +m^^] a = 0. (264) 



To compare with the scalar wave equation we transform the variable as a ^ C = CL\/t-/{t+t) 
We have the following equation for 

C + + ZoeC = 0, (265) 



where 



+ k ( 2i(l + 2A;2m2) _ _ _ A;(- + A^) + yfc^^^ 

It is straightforward to see that when k = the above equation for C becomes exactly the 
same as the scalar wave equation \240{j for k = 0. Therefore in particular the decaying 
properties for the scalar field and odd electromagnetic field are exactly the same in this 
case. 

The even part 

The even part of vector potential can be put in the form 

Aa = aoS{dt)a + aiSa + asSa^a, (267) 
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where aj = ai{t), {I = 0, 1, and 3). It is straightforward to calculate the components of the 
electromagnetic tensor through Eq. (|254() . Noting the basic relations for our eigenfunctions 
OaS = XSa and xsS = —mS, we have '^^ 

FoA = {ai - Xao)SA 

= i^iSA, (268) 
Fo3 = (as + mao)S 

= 1P2S, (269) 
Fa3 = (Aas + mai)SA 

= tp3SA, (270) 
Fab = 0. (271) 

We have defined field strengths of the electric field {ipi and '02) and magnetic field {tp^) 
above, which are all functions of time. Note that they are constrained by definition to 
satisfy 

1P3 = m^pi + Xi^2- (272) 

The strengths ip2, and "03 are invariant under the (mode-decomposed) gauge trans- 
formation 

Aa^ Aa+da{fS), (273) 

where f = f{t) is an arbitrary function of time. It is easy to see this invariance using the 
following induced transformation: 

ao ^ ao + /, 

ai^ai + Xf, (274) 
03 ^ as - mf. 



Next, look at the components of the Maxwell equation (|255)) . which are given by 

(i) /X = 0: 

-A^i" Vi + "1^2 V2 = 0. (275) 

(ii) fi = A: 

^1 = f I - ^ - + - ^Vrni^s. (276) 

\^A^ qi 2q2j Qi 

(iii) /i = 3: 

^2 = f I - ^ - ^^"l ^2 + ^^2 - iV^.r^A^s. (277) 

Note that from the two constraints (|272|) and (|275|) functions ipi and ip2 are solved in 
terms of ^3 as 

V-l = (A'^zr' + m^q^^r'mq^^s, (278) 

^p2 = {X^Qi' + m^q^')-'Xq^^i>3. (279) 



^"As in the odd case, capital Latin indices A, B stand for 1 or 2. 
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Equations (|278|) and (|276|) can virtually be thought of as an unconstrained gauge-invariant 
system of Hamiltonian equations for the variable '03 ^-nd its "momentum" tl^i. If we 
eliminate ■01 we obtain a wave equation 

*^ - - - - I- + 'i^'t>^ ) *3 + iV=(A=,r' + rn'li')*. = 0. (280) 
If using the exact solution we have 

Instead of getting an equation for ip^, we can also get one for ipi, which is 



Since ipi and the "derivative" of ips are related with Ea. H278() . we can reproduce the 
equation for ipi from that of by differentiating Eq. (|281j) (with a further conformal 
transformation). On the other hand, there is no function f{t) such that only the conformal 
transformation ips ^ ip = fip3 reproduces the equation for tpi. We remark that it is the 
equation for ■01 that has a direct resemblance to the odd equation (|264j) , since they coincide 
when k = 0. 
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